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Abstract Starting from the earlier notions of stationary action principles, 
we show how Julian Schwinger’s Quantum Action Principle descended from 
Dirac’s formulation, which independently led Feynman to his path-integral 
formulation of quantum mechanics. The connection between the two is 
brought out, and applications are discussed. The Keldysh-Schwinger time- 
cycle method of extracting matrix elements in nonequilibrium situations is 
described. The variational formulation of quanum field theory and the de¬ 
velopment of source theory constitute the latter part of this work. In this 
document, derived from Schwinger’s lectures over four decades, the continu¬ 
ity of concepts, such as that of Green’s functions, becomes apparent. 
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Chapter 1 

Historical Introduction 


Variational principles for dynamical systems have a long history. Although 
precursors go back at least to Leibnitz (see for example [Euler 175^ 1 and 
Euler [Euler 1744] the “principle of least action” was given modern form by 
de Maupertuis [Maupertuis 1744) [Maupertuis 1746] . We will not attempt to 
trace the history here; a brief useful account is given in Sommerfeld’s lectures 
[Sommerfeld 1964] . The most important names in the history of the develop¬ 
ment of dynamical systems, or at least those that will bear most directly on 
the following discussion, are those of Joseph-Lo uis Lagrange [Lagrange 1788| 
and William Rowan Hamilton [Hamilton 18341 [Hamilton 1835| . 

Here we are concentrating on the work of Julian Schwinger (1918-1994), 
who had profound and pervasive influence on 20th century physics, and whose 
many students have become leaders in diverse fields0 For biographical infor¬ 
mation about his life and work see [Mehra 20001 IMilton 2007] . Therefore, 
we will take up the story in the modern era. Shortly after Dirac’s work 
with Fock and Podolsky [Dirac 1932] . in which the demonstration of the 
equivalence between his theory of quantum electrodynamics, and that of 
Heisenberg and Pauli [Heisenberg 1929| , P. A. M. Dirac wrote a paper on 
“The Lagrangian in Quantum Mechanics” [Dirac 19^ . This paper had a 
profound influence on Richard Feynman’s doctoral dissertation at Princeton 
on “The Principles of Least Action in Quantum Mechanics” [Feynman 1942| , 
and on his later work on the formulations of the “Space-Time Approach to 
Quantum Electrodynamics” [Feynman 1949| . Dirac’s paper further formed 
the basis for Schwinger’s development of the quantum action principle, which 
first appeared in his final operator field formulation of quantum field theory 
[Schwinger 1951b , 


which we will describe in Chapter 6. 


The response of Feynman and Schwinger to Dirac’s inspiring paper was 
completely different. Feynman was to give a global “solution” to the problem 
of determining the transformation function, the probability amplitude con- 


^ For complex reasons, Schwinger’s influence on modern physics is not widely appre¬ 
ciated. His contributions to our current understanding of nature are underrepresented 
in textbooks, with some notable exceptions IToms 20071 . 
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1 Historical Introduction 


necting the state of the system at one time to that at a later time, in terms 
of a sum over classical trajectories, the famous path integral. Schwinger, 
instead, derived (initially postulated) a differential equation for that trans¬ 
formation function in terms of a quantum action functional. This differential 
equation possessed Feynman’s path integral as a formal solution, which re¬ 
mained poorly defined; but Schwinger believed throughout his life that his 
approach was “more general, more elegant, more useful, and more tied to the 
historical line of development as the quantum transcription of Hamilton’s 
action principle” [Schwinger 1973| . 

Later, in a tribute to Feynman, Schwinger commented further. Dirac, of 
course, was the father of transformation theory [Dirac 1927] . The transfor¬ 
mation function from a description of the system at time t 2 to a description 
at time ti is “the product of all the transformations functions associated 
with the successive infinitesimal increments in time.” Dirac said the latter, 
that is, the transformation function from time t to time t + dt corresponds 
to exp[{i/h)dt L], where L is the Lagrangian expressed in terms of the co¬ 
ordinates at the two times. For the transformation function between t 2 and 
ti “the integrand is exp[(i/?i)VF]. where W = dtL.” “Now we know, and 
Dirac surely knew, that to within a constant factor the ‘correspondence,’ 
for infinitesimal dt, is an equality when we deal with a system of nonrela- 
tivistic particles possessing a coordinate-dependent potential energy V .... 
Why then, did Dirac not make a more precise, if less general statement? Be¬ 
cause he was interested in a general question: What, in quantum mechanics, 
corresponds to the classical principle of stationary action?” 

“Why, in the decade that followed, didn’t someone pick up the compu¬ 
tational possibilities offered by this integral approach to the time transfor¬ 
mation function? To answer this question bluntly, perhaps no one needed 
it—until Feynman came along.” [Schwinger 1989a1 . 

But Schwinger followed the differential route, and starting in early 1950 
began a new, his third, formulation of quantum electrodynamics, based on 
a variational approach. This was first published in 1951 [Schwinger 1951b . 
A bit later he started developing a new formulation of quantum kinematics, 
which he called Measurement Algebra, which got its first public presentation 
at Ecole de Physique des Houches in the summer of 1955. There were several 
short notes in the Proceedings of the US National Academy published in 
1960, explaining both the quantum kinematical approach and the dynam¬ 
ical action principle [Schwinger lOGOa] [Schwinger 1960b] [Schwinger 1960c 
[Schwinger 1960d| , but although he often promised to write a book on the 
subject (as he also promised a book on quantum held theory) nothing came 
of it. Les Houches lectures, based on notes taken by Robert Kohler, eventually 
appeared in 1970 [Schwinger 1970a[ . Lectures based on a UCLA course on 
quantum mechanics by Schwinger were eventually published under Englert’s 
editorship [Schwinger 2001[ . The incompleteness of the written record may 
be partly alleviated by the present essay. 

We start on a classical footing. 
























Chapter 2 

Review of Classical Action Principles 


This section grew out of lectures given by Schwinger at UCLA around 1974, 
which were substantially transformed into Chapter 8 of Classical Electrody¬ 
namics [Schwinger 1998| . (Remarkably, considering his work on waveguide 
theory during World War II, now partially recorded in Ref. [Milton 2006j , he 
never gave lectures on this subject at Harvard after 1947.) 

We start by reviewing and generalizing the Lagrange-Hamilton principle 
for a single particle. The action, IU 12 , is defined as the time integral of the 
Lagrangian, L, where the integration extends from an initial configuration or 
state at time ^2 to a final state at time t \: 

Wi2=f"dtL. ( 2 . 1 ) 

Jt2 

The integral refers to any path, any line of time development, from the initial 
to the final state, as shown in Fig. 12. II The actual time evolution of the system 



Fig. 2.1 A possible path from initial state to final state. 
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2 Review of Classical Action Principles 


is selected by the principle of stationary action: In response to infinitesimal 
variations of the integration path, the action W 12 is stationary—does not 
have a corresponding infinitesimal change—for variations about the correct 
path, provided the initial and final configurations are held fixed. 


SW 12 = 0. (2.2) 

This means that, if we allow infinitesimal changes at the initial and final 
times, including alterations of those times, the only contribution to SW 12 
then comes from the endpoint variations, or 


SWi2 = Gi-G2, (2.3) 

where Ga, a = 1 or 2, is a function, called the generator, depending on 
dynamical variables only at time ta- In the following, we will consider three 
different realizations of the action principle, where, for simplicity, we will 
restrict our attention to a single particle. 


2.1 Lagrangian Viewpoint 

The nonrelativistic motion of a particle of mass m moving in a potential 
V (r, t) is described by the Lagrangian 

1 /dr V 

Here, the independent variables are r and t, so that two kinds of variations 
can be considered. First, a particular motion is altered infinitesimally, that 
is, the path is changed by an amount dr: 

r(t) — r(t) + Sr(t). (2.5) 

Second, the final and initial times can be altered infinitesimally, by 6ti and 
St 2 , respectively. It is more convenient, however, to think of these time dis¬ 
placements as produced by a continuous variation of the time parameter, 

t —^ t -j- (5t(t), (2.6) 

so chosen that, at the endpoints. 


St{ti) = 6ti, St{t2) = St2. 
The corresponding change in the time differential is 


(2.7) 


2.1 Lagrangian Viewpoint 


dSt 


dt —^ d(t St^ — j dt^ 

which implies the transformation of the time derivative, 

d / dSt\ d 
dt V dt ) dt 


( 2 . 8 ) 


(2.9) 


Because of this redefinition of the time variable, the limits of integration in 
the action. 


Wi2 = 


1 (drf , 

-m- - dtV 

2 dt 


( 2 . 10 ) 


are not changed, the time displacement being produced through St(t) subject 
to (12.71) . The resulting variation in the action is now 


SWi2 = f dtl • -^Sr - Sr ■ VV - ^ 
.2 1 dt dt dt 


-m 


■V 




= L * Tt 


dr 


1 


m— "dr — -m — +V ] St 


dt 


dr 


dt 


Sr • 


—TO—rr — W 

dt^ 


St 


dt 


V 


-V 

dt , 


( 2 . 11 ) 


where, in the last form, we have integrated by parts in order to isolate dr and 
St. 

Because dr and St are independent variations, the principle of stationary 
action implies that the actual motion is governed by 




d 

dt 


'1 

UrV 

-TO 


2 

KdtJ \ 


= ^V 
dt ’ 


while the total time derivative gives the change at the endpoints, 

G = p • dr — ESt, 

with 


(2.12a) 

(2.12b) 

(2.12c) 


dr 


1 


momentum = p = energy = E = —to I -p ) +V. (2.12d) 


dt 


dr 


dt 


Therefore, we have derived Newton’s second law [the equation of motion 
in second-order form], (l2.12aL and, for a static potential, dV/dt = 0, the 
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2 Review of Classical Action Principles 


conservation of energy, (I2.12bl) . The significance of (I2.12cp will be discussed 
later in Section [231 


2.2 Hamiltonian Viewpoint 


Using the above definition of the momentnm, we can rewrite the Lagrangian 
as 

dv 

L = p- — - H{r,p,t), (2.13) 

where we have introduced the Hamiltonian 


We are here to regard r, p, and t as independent variables in 


Wi 2 = j [p-dr — dtH], 


(2.14) 


(2.15) 


The change in the action, when r, p, and t are all varied, is 


SWi2 = J dt 
= f dt 


d dH dr dH d5t dH 


dt 


+(5p ■ 


dr dH 
dt dp 


St 


dt dr 
dH dH 
dt dt 


The action principle then implies 


= -VH, 


dr dH p 

dt dp TO ’ 

dp dH 

dt dr 

dH _ dH 
dt dt ’ 

G = p ■ (5r — HSt. 


(2.16) 

(2.17a) 

(2.17b) 

(2.17c) 

(2.17d) 


In contrast with the Lagrangian differential equations of motion, which 
involve second derivatives, these Hamiltonian equations contain only first 
derivatives; they are called first-order equations. They describe the same 
physical system, because when (I2.17al) is substituted into (I2.17b|) . we re¬ 
cover the Lagrangian-Newtonian equation (I2.12a|) . Furthermore, if we insert 
(|2.17al) into the Hamiltonian (I2.14L we identify H with E. The third equa- 




















2.3 A Third, Schwingerian, Viewpoint 
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tion (|2.17cl) is then identical with (I2.12bl) . We also note the equivalence of 
the two versions of G. 

But probably the most direct way of seeing that the same physical system 
is involved comes by writing the Lagrangian in the Hamiltonian viewpoint as 


L 


m 

T 


dr\ 

dt J 


-V- 


2m 


p — TO 


dr 

dt 


(2.18) 


The result of varying p in the stationary action principle is to produce 

dv 

p = TO—. (2.19) 

But, if we accept this as the definition of p, the corresponding term in L dis¬ 
appears and we explicitly regain the Lagrangian description. We are justified 
in completely omitting the last term on the right side of (I2.18|) . despite its 
dependence on the variables r and t, because of its quadratic structure. Its 
explicit contribution to SL is 


1 

TO 




- d ^ dr dSt \ 

Jp-m-Jr + m-—J, 


( 2 . 20 ) 


and the equation supplied by the stationary action principle for p variations, 
(|2.19|) . also guarantees that there is no contribution here to the results of r 
and t variations. 


2.3 A Third, Schwingerian, Viewpoint 

Here we take r, p, and the velocity, v, as independent variables, so that the 
Lagrangian is written in the form 

/ fjj' \ 1 

L = p- y— - Vj -b -mv^ -V{r,t) = p- — - H{r, p, v, t), (2.21) 

where 

iJ(r, p, v,t) = p • V - + V{r,t). (2.22) 

The variation of the action is now 
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2 Review of Classical Action Principles 


SW] 


12 


[p • dr — H dt] 


dt 


\ dr d 


dH ^ dH ^ dH 

— -(IP- ^--5“ 

or op 0 v 


_ 

dt dt 


dt 


d, j j [dp dH'' 


dH 


so that the action principle implies 


dp dH „ 

-f- =-= -vy 

dt dr 

dr dH 
dt dp ’ 
dH 

0 =-— = -P + TOV, 

dv 

dH _ dH 
dt dt 
G = p • 6r — HSt. 


(2.24a) 

(2.24b) 

(2.24c) 

(2.24d) 

(2.24e) 


Notice that there is no equation of motion for v since dv/dt does not occur 
in the Lagrangian, nor is it multiplied by a time derivative. Consequently, 
(|2.24cl) refers to a single time and is an equation of constraint. 

From this third approach, we have the option of returning to either of 
the other two viewpoints by imposing an appropriate restriction. Thus, if we 
write (|2.22l) as 


■f^(r,P,v,t) 


+ C(r,t) - ^(p - mv)2, 


(2.25) 


and we adopt 


V = —p 

m 


(2.26) 


as the definition of v, we recover the Hamiltonian description, (12.131) and 
p.l4|) . Alternatively, we can present the Lagrangian ()2.2ip as 


dr 


dt 


^ = ^ 37 -^^ + (P-"iv). —-V -- —-V 


dr 


dt 


dr 


dt 


(2.27) 


Then, if we adopt the following as definitions, 
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dr 

v=-, p = mv, 


(2.28) 


the resultant form of L is that of the Lagrangian viewpoint, (12.411 . It might 
seem that only the definition v = dr/dt, inserted in (I2.27p . suffices to regain 
the Lagrangian description. But then the next to last term in (12.271) would 
give the following additional contribution to SL, associated with the variation 
Sr: 


(p — mv) • —(5r. 
' dt 


(2.29) 


In the next Chapter, where the action formulation of electrodynamics is 
considered, we will see the advantage of adopting this third approach, which 
is characterized by the introduction of additional variables, similar to v, for 
which there are no equations of motion. 


2.4 Invariance and Conservation Laws 

There is more content to the principle of stationary action than equations of 
motion. Suppose one considers a variation such that 


6Wi2 = 0, (2.30) 

independently of the choice of initial and final times. We say that the action, 
which is left unchanged, is invariant under this alteration of path. Then the 
stationary action principle (j2.3|) asserts that 

SWi2 = Gi - G2 = 0, (2.31) 

or, there is a quantity G(t) that has the same value for any choice of time t; 
it is conserved in time. A differential statement of that is 

= 0. (2.32) 

The G functions, which are usually referred to as generators, express the 
interrelation between conservation laws and invariances of the system. 

Invariance implies conservation, and vice versa. A more precise statement 
is the following: 

If there is a conservation law, the action is stationary under an infinitesimal 
transformation in an appropriate variable. 

The converse of this statement is also true. 

If the action W is invariant under an infinitesimal transformation (that is, 
5W = 0), then there is a corresponding conservation law. 


This is the celebrated theorem of Amalie Emmy Noether [Noether 1918] . 
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2 Review of Classical Action Principles 


Here are some examples. Suppose the Hamiltonian of (12.131) does not de¬ 
pend explicitly on time, or 

Wi 2 = J [p- dr - H{r,p)dt]. (2.33) 

Then the variation (which as a rigid displacement in time, amounts to a shift 
in the time origin) 

6t = constant (2.34) 

will give 5Wi2 = 0 [see the first line of (I2.16p . with i5r = 0, (5p = 0, dSt/dt = 0, 
dH/dt = 0]. The conclusion is that G in (|2.17dl) . which here is just 

Gt = -H6t, (2.35) 


is a conserved quantity, or that 



(2.36) 


This inference, that the Hamiltonian—the energy—is conserved, if there is 
no explicit time dependence in H, is already present in (I2.17c[) . But now a 
more general principle is at work. 

Next, consider an infinitesimal, rigid rotation, one that maintains the 
lengths and scalar products of all vectors. Written explicitly for the posi¬ 
tion vector r, it is 

i5r = (5u;xr, (2.37) 

where the constant vector 6u: gives the direction and magnitude of the rota¬ 
tion (see Fig. 12.21) . Now specialize (I2.14|) to 



Fig. 2.2 Slj X r is perpendicular to Suj and r, and represents an infinitesimal rotation 
of r about the Soj axis. 
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H = ^ + V(r), 


where r = |r|, a rotationally invariant structure. Then 

14^12 = J [p‘ dr — H dt] 


(2.38) 


(2.39) 


is also invariant under the rigid rotation, implying the conservation of 

= P • • rxp. ( 2 . 40 ) 

This is the conservation of angular momentum, 

L = rxp, ^L = 0 . ( 2 . 41 ) 

Of course, this is also contained within the equation of motion, 

4-L =-rx Vy =-rxf^ = 0, (2.42) 

dt or 

since V depends only on |r|. 

Conservation of linear momentum appears analogously when there is in¬ 
variance under a rigid translation. For a single particle, (I2.17bll tells us im¬ 
mediately that p is conserved if F is a constant, say zero. Then, indeed, the 
action 


Wi2 = 


/2 L 


p • dr — ^—dt 
2m 


is invariant under the displacement 

Sr = Se = constant, 

and 

Gse = p • <5e 


(2.43) 


( 2 . 44 ) 


( 2 . 45 ) 


is conserved. But the general principle acts just as easily for, say, a system 
of two particles, a and b, with Hamiltonian 


H=^ + ^ + V{r, - r,). 
Zma 2mh 

This Hamiltonian and the associated action 


(2.46) 


Wi2 = J \Pa' dra + Pb ■ drb - H dt] (2.47) 


are invariant under the rigid translation 
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2 Review of Classical Action Principles 


= Sri, = Se, 


(2.48) 


with the implication that 

Gse = Pa • Sra + p& • Sri, = (pa + Pf,) • Se (2.49) 

is conserved. This is the conservation of the total linear momentum, 

P = Pa + P&, = 0- (2-50) 

at 

Something a bit more general appears when we consider a rigid translation 
that grows linearly in time: 


Sra = Sri, = Svt, (2-51) 

using the example of two particles. This gives each particle the common 
additional velocity Sv, and therefore must also change their momenta, 

(5pa = maSv, Spb = mbSv. (2.52) 

The response of the action (12.471) to this variation is 

SWi2 = J [(Pq + Pb) 'S^dt + Sv (madra + rubdrb) - (Po + Pb) • Sv dt] 

1 

d[(mara+mbrb)’Sv\. (2.53) 

The action is not invariant; its variation has end-point contributions. But 
there is still a conservation law, not of G = P ■ Svt, but of N ■ Sv, where 



N = Pt - (ruara 4- mbrb). 


Written in terms of the center-of-mass position vector 


R = 


rUara + mbrb 

M 


M = ma + mb. 


the statement of conservation of 


(2.54) 


(2.55) 


N = Pt - MR, 


(2.56) 


namely 


0 = 



M- 


.dR 
dt ’ 


(2.57) 


is the familiar fact that the center of mass of an isolated system moves at the 
constant velocity given by the ratio of the total momentum to the total mass 
of that system. 
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2.5 Nonconservation Laws. The Virial Theorem 


The action principle also supplies useful nonconservation laws. Consider, for 
constant i5A, 

(5r = (5Ar, (5p = —i5Ap, (2.58) 

which leaves p • dr invariant, 

(5(p • dr) = (—(5Ap) ■ dr + p • (SXdr) = 0. (2.59) 

But the response of the Hamiltonian 

H = Tip) + Vir), T{p) = ^, (2.60) 


is given by the noninvariant form 

SH = SX{-2T + r-VV). (2.61) 

Therefore we have, for an arbitrary time interval, for the variation of the 
action (12.151) . 


dWi 2 = [ dt[SX{2T-r-W)] = Gi-G 2 = [ dt^(p-dAr) (2.62) 
J 2 J 2 dt 

or, the theorem 

4-r-p = 2T-r-VL. (2.63) 

dt 

For the particular situation of the Coulomb potential between charges, 
V = constant/r, where 


r- VC = r-^V = —C, (2.64) 

dr 

the virial theorem asserts that 

^(r.p) = 2T + C. (2.65) 

We apply this to a bound system produced by a force of attraction. On 
taking the time average of (12.651) the time derivative term disappears. That 
is because, over an arbitrarily long time interval r = L — ^ 2 , the value of 
r • p(ti) can differ by only a finite amount from r • p(t 2 )j and 




fn dt T 


( 2 . 66 ) 


as r —>■ 00 . The conclusion, for time averages. 
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2T = -V, 


(2.67) 


is familiar in elementary discussions of motion in a l/r potential. 

Here is one more example of a nonconservation law: Consider the variations 


Sr = S\- 


Sp = -SA{P-^-^ 


= sx 


rx (rxp) 


Again p • dr is invariant: 
(5(p ■ dr) = —SX f — — 


, , dr ^, r ■ dr 

dr + p ■ ( (5A-(5Ar —— 


and the change of the Hamiltonian (12.601) is now 


SH = SX 


r 

-• VH 


mj-^ r 


The resulting theorem, for H = H(r), is 


sO-) 


dr ’ 


(2.68a) 

(2.68b) 

0, (2.69) 

(2.70) 

(2.71) 


which, when applied to the Coulomb potential, gives the bound-state time 
average relation 



(2.72) 


This relation is significant in hydrogen fine-structure calculations (for exam¬ 
ple, see [Schwinger 2001| ). 
















Chapter 3 

Classical Field Theory—Electrodynamics 


This Chapter again grew out of Schwinger’s UCLA lectures. These evolved, 
torturously, into Chapter 9 of Classical Electrodynamics [Schwinger 1998 . 
Here we use Gaussian units. 


3.1 Action of Particle in Field 


It was stated in our review of mechanical action principles in the previous 
Chapter that the third viewpoint, which employs the variables r, p, and v, 
was particularly convenient for describing electromagnetic forces on charged 
particles. With the explicit, and linear, appearance of v in what plays the 
role of the potential function when magnetic fields are present, we begin to 
see the basis for that remark. Indeed, we have only to consult (I2.2ip to find 
the appropriate Lagrangian: 

where (j) and A are the scalar and vector potentials, respectively. To recapit¬ 
ulate, the equations resulting from variations of p, r, and v are, respectively. 


dr 

d „ 
TtO = 


'-V • A 

c 


p = mv H—A. 

c 


(3.2a) 

(3.2b) 

(3.2c) 


We can now move to either the Lagrangian or the Hamiltonian formulation. 
For the first, we simply adopt v = dr/dt as a definition (but see the discussion 
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3 Classical Field Theory—Electrodynamics 


in Sec. lOl) and get 



• A. 


(3.3) 


Alternatively, we use (I3.2cp to define 


V 



(3.4) 


and find 



(3.5b) 


(3.5a) 


3.2 Electrodynamic Action 

The electromagnetic field is a mechanical system. It contributes its variables 
to the action, to the Lagrangian of the whole system of charges and fields. 
In contrast with the point charges, the field is distributed in space. Its La¬ 
grangian should therefore be, not a summation over discrete points, but an 
integration over all spatial volume elements. 



(3.6) 


this introduces the Lagrange function, or Lagrangian density, £. The total 
Lagrangian must be the sum of the particle part, dm, and the field part, 
(I3lil) . where the latter must be chosen so as to give the Maxwell equations, 
in Gaussian units: 



(3.7a) 


V ■ E = 47rp, 


(3.7b) 


V-B = 0. 


The homogeneous equations here are equivalent to the construction of the 
electromagnetic field in term of potentials, or, 



E-Vc/>, 


(3.8a) 

(3.8b) 


B = VxA. 




3.2 Electrodynamic Action 
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Thus, we recognize that A(r,t), E(r,t), in analogy with r(t), p(t), obey 
equations of motion while (j){v, t), B(r, t), as analogues of v(t), do not. There 
are enough clues here to give the structure of £fieidi apart from an overall 
factor. The anticipated complete Lagrangian for microscopic electrodynamics 
is 


^ = E 

a 

1 , 


p« • ( ^ - '"a ) + ^^avl - ea(t){ra) + —• A(r^) 


+ 




(dr) 


dt 


E- - 


l_a 

'.dt^ 


A-V0 -B-VxA + -(b2-£;2) 


.(3.9) 


The terms that are summed in (EH) describe the behavior of charged 
particles under the influence of the fields, while the terms that are integrated 
describe the field behavior. The independent variables are 


ra(t), Va(t), Pa(t), '('(r,t), A(r,t), E(r,t), 


B(r,t), t. (3.10) 


We now look at the response of the Lagrangian to variations in each of these 
variables separately, starting with the particle part: 


(5ro : 5L = —{Sra-Pa) + 5va 


dVa 

dt 




dVa : 5L = (5vq 
(5pq : 5L = dpa 


-Pa + rnaVa + a) 

dVg _ 
dt 



(3.11a) 

(3.11b) 

(3.11c) 


The stationary action principle now implies the equations of motion 


-^ =-eaVa ((/>(ra) - ^ • A(r^)) , (3.12a) 

m-aVo = Pa - ^A(r^), (3.12b) 

Va = (3.12c) 


which are the known results, (I3.2all ~ ()3.2c|) . 

The real work now lies in deriving the equations of motion for the fields. 
In order to cast all the field-dependent terms into integral form, we introduce 
charge and current densities. 


Pir,t) = '^ead{r-ra{t)), (3.13a) 

a 

j(j', 0 = E eaVa(t)(5(r - rg{t)), 


(3.13b) 
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SO that 




-ea(l){ra) + -^Va ■ A(r^) = / (dr) 


+ -j • A 

c 


(3.14) 


The volume integrals extend over sufficiently large regions to contain all the 
fields of interest. Consequently, we can integrate by parts and ignore the 
surface terms. The responses of the Lagrangian (j3.9|) to field variations, and 
the corresponding equations of motion deduced from the action principle are 


d(/) : 

SL 


V-E 

dA : 

SL 


VxB 

dE : 

SL 


E 

dB : 

SL 


SL = — J (dr) i5^(V • E — 47rp), 


1 d 
47rc dt 

1 d „ 


(dr) dA ■ 
47r. 


1 (9E 47r. „ „ 

—J-VxB 

c ot c 


(3.15a) 
(3.15b) 

(3.15c) 

(3.15d) 

(3.15e) 

(3.15f) 

(3.15g) 
(3.15h) 

We therefore recover Maxwell’s equations, two of which are implicit in the 
construction of E and B in terms of potentials. By making a time variation 
of the action [variations due to the time dependence of the fields vanish by 
virtue of the stationary action principle—that is, they are already subsumed 
in Eqs. (13.151) ]. 


-i^A-V(^-E 

c ot 


cdt c 

^ I {dr) St 

-~A-V0, 

c at 

SL= — f (dr) dB • (-Vx A + B), 
47r J 

B = VxA. 




St : 5W = J dt 

we identify the Hamiltonian of the system to be 

H = Y1 (P“ “ • W - ^rUavl + ea</>(ra) 


(3.16) 


a ^ 
1 


/ <*> 


E • V(/> + B • Vx A + -{E^ - B'^) 


(3.17) 
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which is a constant of the motion, dH/dt = 0. The generators are inferred 
from the total time derivative terms in (|3.11ap . (I3.15cl) . and (13.161) . 


to be 


5W\2 — Gi — G2, 


(3.18a) 


G = Jrg • Pa 
a 


1 

47rc 



E ■ ,5A - HSt. 


(3.18b) 


3.3 Energy 


Notice that the total Lagrangian (13.9p can be presented as 

i = (3.19) 


where the Hamiltonian is given by (I3.17|) . The narrower, Hamiltonian, de¬ 
scription is reached by eliminating all variables that do not obey equations of 
motion, and, correspondingly, do not appear in G. Those “superfluous” vari¬ 
ables are the Va and the fields (j) and B, which are eliminated by using (I3.12bp , 
(j3.15bll . and (I3.15hll . the equations without time derivatives, resulting, first, 
in the intermediate form 


^ = E 

a 






'E^ + B^ 
Stt 



(3.20) 

The first term here is the energy of the particles moving in the field [par¬ 
ticle energy—see (j3.5bl) ]. so we might call the second term the field energy. 
The ambiguity of these terms (whether the potential energy of particles is 
attributed to them or to the fields, or to both) is evident from the existence 
of a simpler form of the Hamiltonian 


H = 





E^ + B^ 
87r 


B = V X A, (3.21) 


where we have used the equivalence of the two terms involving (j), given in 

(EH. 

This apparently startling result suggests that the scalar potential has dis¬ 
appeared from the dynamical description. But, in fact, it has not. If we vary 
the Lagrangian (j3.19() , where H is given by (j3.21() , with respect to E we find 

3i=_^/(*)3E.(i|A + E)=0. 


(3.22) 
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Do we conclude that + E = 0? That would be true if the <5E(r, t) were 

arbitrary. They are not; E is subject to the restriction—the constraint— 
(|3.15bl) . which means that any change in E must obey 


V-5E = 0. 


(3.23) 


The proper conclusion is that the vector multiplying (5E in (13.221) is the 
gradient of a scalar function, just as in (I3.15fj) . 

i^A + E = -V</), (3.24) 

c ot 

for that leads to 

SL = -^ J{dr){V = (3.25) 

as required. 

The fact that the energy is conserved, 



(3.26) 


where 

H = ^^rnavl +J{dr)U, U , (3.27) 

is a simple sum of particle kinetic energy and integrated field energy density, 
can be verified directly by taking the time derivative of (13.201) . The time rate 
of change of the particle energy is computed directly: 

^ + ea(/>(ra)^ • (3-28) 

a ^ ''a 


We can compute the time derivative of the field energy by using the equation 
of energy conservation, 

jJidr)U = -Jidr)3‘B, (3.29) 

to be 



(3.30) 
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Here we have used (I3.15f|) . and have noted that 


(dr) 




= 0 


(3.31) 


by charge conservation. Observe that (13.281) and (I3.30|) are equal in magnitude 
and opposite in sign, so that their sum is zero. This proves the statement of 
energy conservation (13.26(1 . 


3.4 Momentum and Angular Momentum Conservation 

The action principle not only provides us with the field equations, particle 
equations of motion, and expressions for the energy, but also with the gen¬ 
erators (I3.18bl) . The generators provide a connection between conservation 
laws and invariances of the action (recall Section . Here we will further 
illustrate this connection by deriving momentum and angular momentum 
conservation from the invariance of the action under rigid coordinate trans¬ 
lations and rotations, respectively. [In a similar way we could derive energy 
conservation, (13.261) . from the invariance under time displacements—see also 
Section 13.6| . 

Under an infinitesimal rigid coordinate displacement, de, a given point 
which is described by r in the old coordinate system is described by r -|- 5e 
in the new one. (See Fig. 13.11 1 The response of the particle term in (I3.18b|) 


new old 



Fig. 3.1 Rigid coordinate displacement, where the new coordinate system is dis¬ 
placed by a rigid translation —relative to the old coordinate system. 


is simple: (5e • X^aPai field part, we require the change, 5A., of the 

vector potential induced by the rigid coordinate displacement. The value of 
a field at a physical point P is unchanged under such a displacement, so 
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that if r and r + are the coordinates of P in the two frames, there are 
corresponding functions F and F such that 

F{P) = F{r)=F(r + Se), (3.32) 

that is, the new function F of the new coordinate equals the old function F 
of the old coordinate. The change in the function F at the same coordinate 
is given by 

F{v) = F{v) + 5F{v), (3.33) 

so that 

5F{v) = F{y — Se) — F(r) = —Se ■ VF{r), (3.34) 

for a rigid translation (not a rotation). 

As an example, consider the charge density 

p(r) = ^ea(5(r - ra). (3.35) 

a 

If the positions of all the particles, the ra, are displaced by i5e, the charge 
density changes to 


p(r) + (5p(r) = ^ea(5(r - Ta - (5e), (3.36) 

a 

where 

(5(r - Ta - Se) = (5(r - r^) - Se • Vr(5(r - r^), (3.37) 

and therefore 

in agreement with p.34l) . 

So the field part of G in 

-/(*)jEE.iA = 


(3.39) 

where the last rearrangement makes use of (I3.15b() and (j3.15hl) , and the vector 
identity 

<5ex(VxA) = V(^e-A)-((5e-V)A. (3.40) 

Including the particle part from (|3.18bll we find the generator corresponding 
to a rigid coordinate displacement can be written as 


Sp{r) = -Se-Vp{r), 


(3.38) 


(I3.18b() is 

^ J{dr)E,{Se-V)A, 
i eaSe • A(r^) + ^ f (dr) (E X B) • de. 


G = (5e • P 


(3.41) 
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where 

P = ^ (pa - ^A(r^)) f (dr) ExB = ^ ma^^a+ f {dr) G, (3.42) 

r, n 


with G the momentum density. Since the action is invariant under a rigid 
displacement, 

0 = dtE = Gi-G2 = (Pi-P2)-<5r, (3.43) 

we see that 

Pi = P2, (3.44) 

that is, the total momentum, P, is conserved. This, of course, can also be 
verified by explicit calculation: 


_d 

dt 


(dr)—ExB 



pE + - j X B 

c 


+ ^VaXB(rJ^ , 

n. ^ ' 


(3.45) 


from which the constancy of P follows. 

Similar arguments can be carried out for a rigid rotation for which the 
change in the coordinate vector is 


dr = duixr, 


(3.46) 


with biiS constant. The corresponding change in a vector function is 

A(r + dr) = A(r) + da;xA(r) (3.47) 


since a vector transforms in the same way as r, so the new function at the 
initial numerical values of the coordinates is 

A(r) = A(r) — (dr • V)A(r) + dwx A(r). (3.48) 

The change in the vector potential is 

dA = -(dr. V)A + da;xA. (3.49) 

The generator can now be written in the form 

G = dw • J, (3.50) 


where the total angular momentum, J, is found to be 
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which again is a constant of the motion. 


3.5 Gauge Invariance and the Conservation of Charge 


An electromagnetic system possesses a conservation law, that of electric 
charge, which has no place in the usual mechanical framework. It is con¬ 
nected to a further invariance of the electromagnetic fields—the potentials 
are not uniquely defined in that if we let 


A +VA, 



(3.52) 


the electric and magnetic fields defined by (I3.15fp and (I3.15hl) remain unal¬ 
tered, for an arbitrary function A. This is called gauge invariance; the cor¬ 
responding substitution (I3.52|) is a gauge transformation. [The term has its 
origin in a now obsolete theory of Hermann Weyl (1885-1955) |Weyl 191^ .] 
This invariance of the action must imply a corresponding conservation law. 
To determine what is conserved, we compute the change in the Lagrangian, 
(EH), explicitly. Trivially, the field part of L remains unchanged. In consid¬ 
ering the change of the particle part, we recognize that (13.521) is incomplete; 
since v is a physical quantity, p — (e/c)A must be invariant under a gauge 
transformation, which will only be true if (13.521) is supplemented by 

p—;>p-l--VA. (3.53) 

c 

Under the transformation (13.521) and (13.531) , the Lagrangian becomes 


L L = L -I- 


-VA- 




dVa 

dt 

dt 


ea d ea ^ ' 

—-^AH-Va-VA 

c ot c 




where 


= y]—A(r„,t). 


(3.54) 

(3.55) 


What is the physical consequence of adding a total time derivative to a La¬ 
grangian? It does not change the equations of motion, so the system is unal¬ 
tered. Since the entire change is in the end point behavior. 


ITi 2 = Wi2 -l- {w\ — W2), 


(3.56) 
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the whole effect is a redefinition of the generators, G, 

G = G + 8w. 


(3.57) 


This alteration reflects the fact that the Lagrangian itself is ambiguous up 
to a total time derivative term. [This term may also be seen as arising from 
the field term in the generator (I3.18bl) .] 

To ascertain the implication of gauge invariance, we rewrite the change in 
the Lagrangian given in the first line of (j3.54l) by use of (I3.12cl) . 


L-L 



P^A+j.VA 


(3.58) 


and apply this result to an inhnitesimal gauge transformation, A —>■ (5A. The 
change in the action is then 

SWi 2 = Gsm - Gsx, - dt y (dr) -JX + V • (3.59) 

with the generator being 

Gsx = J{dr)^pSX. (3.60) 


In view of the arbitrary nature of i5A(r, t), the stationary action principle now 
demands that, at every point, 

^p + V.j = 0, (3.61) 

that is, gauge invariance implies local charge conservation. (Of course, this 
same result follows from Maxwell’s equations.) Then, the special situation 
SX = constant, where SA = 64> = 0, and W 12 is certainly invariant, implies a 
conservation law, that of 

G5a = -<5AQ, (3.62) 

c 

in which 

Q = J(dr) p (3.63) 

is the conserved total charge. 


3.6 Gauge Invariance and Local Conservation Laws 

We have just derived the local conservation law of electric charge. Electric 
charge is a property carried only by the particles, not by the electromag- 
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netic field. In contrast, the mechanical properties of energy, linear momen¬ 
tum, and angular momentum are attributes of both particles and fields. For 
these we have conservation laws of total quantities. What about local con¬ 
servation laws? The usual development of electrodynamics refers to local 
non-conservation laws; they concentrated on the fields and characterized the 
charged particles as sources (or sinks) of field mechanical properties. It is 
natural to ask for a more even-handed treatment of both charges and fields. 
We shall supply it, in the framework of a particular example. The property 
of gauge invariance will be both a valuable guide, and an aid to simplifying 
the calculations. 

The time displacement of a complete physical system identifies its total 
energy. This suggests that time displacement of a part of the system provides 
energetic information about that portion. The ultimate limit of this spatial 
subdivision, a local description, should appear in response to an (infinites¬ 
imal) time displacement that varies arbitrarily in space as well as in time. 


St{r, t). 


Now we need a clue. How do fields, and potentials, respond to such coordi¬ 
nate-dependent displacements? This is where the freedom of gauge transfor¬ 
mations enters: The change of the vector and scalar potentials, by VA(r, t), 
— {l/c){d/dt)X{r,t), respectively, serves as a model for the potentials them¬ 
selves. The advantage here is that the response of the scalar A(r,t) to the 
time displacement can be reasonably taken to be 


(A -I- i5A)(r, t + 6t) = A(r, t), 


(3.64a) 


or 



(3.64b) 


Then we derive 



(3.65b) 


(3.65a) 


which is immediately generalized to 



dt 


(3.66a) 



(3.66b) 


or, equivalently. 


i5A = c(5tE -I- V {(j)c6t) 
c ot 


(3.67a) 

(3.67b) 
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In the latter form we recognize a gauge transformation, produced by the scalar 
(j)cSt, which will not contribute to the changes of field strengths. Accordingly, 
for that calculation we have, effectively, 6A = cStE, 5<j) = 0, leading to 

(5E = — — (cSm = (3.68a) 

c ot ot at 

SB = Vx (cStB) = -St^B -Ex VcSt] (3.68b) 

the last line employs the field equation VxE = —{l/c){dB/dt). 

In the following we adopt a viewpoint in which such homogeneous field 
equations are accepted as consequences of the definition of the fields in terms 
of potentials. That permits the field Lagrange function (13.91) to be simplified: 


>Cfield= ^(^2-^2). (3.69) 

Stt 

Then we can apply the field variation (|3.68bp directly, and get 

d \ d c 

SCiield = —St — Cdeld — ~ 

at AtT at dTT 

= - T(e 2 + - ^ExB . VSt. (3.70) 

ot Stt ot 47r 

Before commenting on these last, not unfamiliar, field structures, we turn to 
the charged particles and put them on a somewhat similar footing in terms 
of a continuous, rather than a discrete, description. 

We therefore present the Lagrangian of the charges in (13.91) in terms of a 
corresponding Lagrange function, 

Lcharges 


— / ■^charges; 


(3.71a) 


where 


and 


C 


charges 




(3.71b) 


£a = S{r - ra{t)) 


1 c 

-maVaitf - ea(/)(ra, t) + -^^o.{t) ’ A(r^, t) 


■ (3.71c) 


the latter adopts the Lagrangian viewpoint, with = dva/dt accepted as a 
definition. Then, the effect of the time displacement on the variables ra{t), 
taken as 


(ra + 6ra){t + St) = ra{t), 

Sra(t) = -St{ra,t)Vait), 


(3.72a) 

(3.72b) 
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implies the velocity variation 


Svait) 


-5t{ra,t)—-Va(t) - Va(t) 


dt 


St - 


VSt 


(3.73) 


the last step exhibits both the explicit and the implicit dependences of 
St{ra, t) on t. In computing the variation of ^(ra, t), for example, we combine 
the potential variation given in (|3.66bl) with the effect of i5ra: 

54>{va{t), t) = ~ ‘ Va<(> = (3.74a) 

and, similarly, 

SA{rait), t) = + (jjcVSt - Stva ■ VaA = -St^A. + (jjcVSt. (3.74b) 

The total effect of these variations on Ca is thus 


SCa = -St-^£a+S{r - ra{t)) (^-ruavl - + CacSj ’ 

(3.75a) 

or 

SCa = -J^iStCa) - S(r - ra{t))Ea + Va • V St^ , (3.75b) 

where we see the kinetic energy of the charged particle, 

Ea = ^mavl- (3.76) 


We have retained the particle symbol d/dt to the last, but now, being firmly 
back in the field, space-time viewpoint, it should be written as djdt, referring 
to all t dependence, with r being held fixed. The union of these various 
contributions to the variation of the total Lagrange function is 

SCtot = -^{StCtot) - Utot-^St - Stot • (3.77) 

at at 

where, from (13.701) and (I3.75bl) . 

Utot = ^{E^+ B^) + ^ 5(r - rait))Ea (3.78a) 


and 

Stot = ^ExB + ^ <5(r - rait))Ea^ra, (3.78b) 

a 

are physically transparent forms for the total energy density and total energy 
flux vector. 
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To focus on what is new in this development, we ignore boundary effects 
in the stationary action principle, by setting the otherwise arbitrary dt{r,t) 
equal to zero at ti and ^ 2 - Then, through partial integration, we conclude 
that ^ 

SWi 2 = dtj (dr) St + V • Stot) = 0, (3.79) 

from which follows the local statement of total energy conservation, 

^17tot + V-Stot = 0,. (3.80) 



Chapter 4 

Quantum Action Principle 


This Chapter, and the following one, are based on lectures given by the author 
in quantum field theory courses at the University of Oklahoma over several 
years, based in turn largely on lectures given by Schwinger at Harvard in the 
late 1960s. 

After the above reminder of classical variational principles, we now turn to 
the dynamics of quantum mechanics. We begin by considering the transfor¬ 
mation function {a', t + dt\b', t). Here \b',t) is a state specified by the values 
b' = {6'} of a complete set of dynamical variables B{t), while |a', t + dt) is a 
state specified by values a' = {o'} of a (different) complete set of dynamical 
variables A{t + dt), defined at a slightly later time[i] We suppose that A and 
B do not possess any explicit time dependence—that is, their definition does 
not depend upon t. Here 


{a',t + dt\ = {a',t\U, (4.1) 

where the infinitesimal time translation operator is related to the generator 
of time translations as follows. 


U = l+iG = l-idtH. (4.2) 

The Hamiltonian iJ is a function of dynamical variables, which we write 
generically as x(0i of t explicitly. Thus 

{a\t + dt\b' ,t) = {a ,t\l — idt H{x{t),t)\b' ,t). (4-3) 

We next translate states and operators to time zero: 

{a',t\ = {a'\U{t), \b',t) = U-\t)\b'), (4.4a) 

x{t) = U~\t)xU{t), (4.4b) 


^ Here Schwinger is using his standard notation, designating eigenvalues by primes. 
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where x = x(0)i etc. Then, 

{a' , t + dt\b', t) = {a'\l — idtH (x, t)\b')^ (4-5) 

or, as a differential equation 

6 dyn{a ,t + dt\b',t) = i{a\Sdyn[-dtH]\b') 

= i{a,t + dt\6dyn[-dtH{x{t),t)]\b',t), (4.6) 

where ^dyn corresponds to changes in initial and final times, 6 t 2 and 6 ti, and 
in the structure of i/, SH. [By reintroducing dt in the state on the left in the 
second line, we make a negligible error of 0 {dt‘^).] 

However, we can also consider kinematical changes. To understand these, 
consider a system defined by coordinates and momenta, {( 7 Q(t)}, {pa{t)}, 
a = 1,..., n, which satisfy the canonical commutation relations, 

[qa(t),pbit)] = iSab, (?i = l) (4.7a) 

[qa{t),qb{t)] = [Pa{t),Pb{t)] = 0. (4.7b) 

A spatial displacement Sqa is induced by 

n 

U=l + iGq, Gq = '^PaSqa. (4.8) 

a—1 

In fact {Sqa is a number, not an operator), 

U~^qaU = qa- ^[qa,Gq] 

= qa- Sqa, (4.9) 

while 

U~^PaU =Pa- ^[Pa,Gq] = Pa- (4.10) 

The (dual) symmetry between position and momentum, 

q^p, p^-q, (4.11) 

gives us the form for the generator of a displacement in p: 

Gp = -^qaSpa- (4.12) 

a 

A kinematic variation in the states is given by the generators 


4i„( i = n-( i = ( \^G, 

4i„| ) = !)-!) = -iG\ ) 


(4.13a) 

(4.13b) 
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SO, for example, under a 5q variation, the transformation function changes by 


5q{a' ,t + dt\b' ,t) = i{a' + dt \^ [pa{t + dt)5qa{t + dt) - pa{t)5qa{t)] \b' ,t). 

(4.14) 

Now the dynamical variables at different times are related by Hamilton’s 
equations. 


dpajt) 

dt 


-[pa{t),H{q{t),p{t),t)] 

i 

dH, , 
oqa 


so 

Paid + dt) -Pa{t) = dt ^^')^ = -dt^{t). 

dt (yQa 

Similarly, the other Hamilton’s equation 

dqg _ 
dt dpa 


(4.15) 

(4.16) 


(4.17) 


implies that 


qa{t + dt) 


, , , dH 

qa{t) = dt—{t). 

dPa 


(4.18) 


From this we deduce first the q variation of the transformation function, 


i5q(a', t + dt\b', t) 

= i(a',t + dt\ ^Pa(t)[(5ga(t + dt) - 5qa{t)] - dt- — 5qa{t) + 0{dt^)\b',t) 
„ oqa 


= i{a', t + dt\Sq 


^Pa(t).[qait + dt) 


Qait)] - dtH{q{t),p{t),t) 


\b',t), 

(4.19) 


where the dot denotes symmetric multiplication of the p and q operators. 
For p variations we have a similar result: 
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Sp{a',t + dt\b',t) 

= -i{a' ,t + dt\'^[qa{t + dt)Spa{t + dt) - qa{t)6pa{t)]\b',t) 

a 

= -i{a' ,t + dtiy^ qa{t)[Spait + dt) - Spa{t)] +dt- — {t)Spait)\b',t) 

„ 9Pa 


= i{a', t + dt\Sp 


'^qait).{pa{t + dt) - Pa{t)) - dt H {q{t) , p{t) ,t) 


\b',t). 

(4.20) 


That is, for q variations 

6g{a', t + dt\b\ t) = i(a', t + dt\Sq [dtLq] \b\ t), (4.21a) 

with the quantum Lagrangian 

Lq = '^Pa-qa- H{q,p,t), (4.21b) 


while for p variations 

6 p{a', t + dt\b', t) = i{a, t + dt\6p [dtLp] \b', t), (4.22a) 

with the quantum Lagrangian 

Lp = - '^qa.pa- H{q,p,t). (4.22b) 

a 

We see here two alternative forms of the quantum Lagrangian. Note that the 
two forms differ by a total time derivative, 

Lq ~ Lp = ^ ^ '^ Pa-qa- (4.23) 

a 

We now can unite the kinematic transformations considered here with the 
dynamic ones considered earlier, in Eq. (14.611 : 

^ = i^dyn + i^kin : d{a ,t + dt\b', dt) = i{a\t + dt\S[dt L]\b',t). (4.24) 


Suppose, for concreteness, that our states are defined by values of g, so 
that 

Sp{a',t + dt\b't) = 0. (4-25) 

This is consistent, as a result of Hamilton’s equations, 


a 


qa - 


dH 

dpa 


= 0 . 


dpLq ^ ( dpf 


(4.26) 
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In the following we will use Lg. 

It is immediately clear that we can iterate the infinitesimal version (I4.24p 
of the quantum action principle by inserting at each time step a complete 
set of intermediate states (to simplify the notation, we ignore their quantum 
numbers): 

(tilh) = - dt){ti - dt\ti - 2dt) • • • (t2 + 2dt\t2 + dt){t2 +dt\t2), (4.27) 


So in this way we deduce the general form of Schwinger’s quantum action 
principle: 

5{ti\t2) = i{ti\5 f dtL\t2)- (4.28) 

Jt2 

This summarizes all the properties of the system. 

Suppose the dynamical system is given, that is, the structure of H does 
not change. Then 

^(^11^2) = *(^1 |Gi — 1^21^2), (4.29) 

where the generator Ga depends on p and q at time ta- Comparing with the 
action principle (14.281) we see 


sf dtL = Gi-G2, (4.30) 

Jt2 


which has exactly the form of the classical action principle (12.31) . except that 
the Lagrangian L and the generators G are now operators. If no changes 
occur at the endpoints, we have the principle of stationary action, 


S 




= 0 . 


(4.31) 


As in the classical case, let us introduce a time parameter t, t = tir), such 
that T 2 and ri are fixed. Calling the new time parameter by the original 
name, the above variation reads 


[5pa-dqa + Pa-dSqa - 6H dt - H dSt] 

a 


= d 


'^PaMa - H 6t 


[5pa-dqa - dpa-Sqa] 


SH dt + dH 5t, 

(4.32) 


so the action principle says 
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G = Y^Pa-ha - H St, 




dt 


dqa 


dpa 


dt 


dt 


(4.33a) 

(4.33b) 


We will again assume Spa, Sqa are not operators (that is, they are proportional 
to the unit operator); then we recover Hamilton’s equations, 


dH _ dH 
dt dt 
dH _ ^ 

dpa dt ’ 

dH _ dpa 
dqa dt 


(4.34a) 

(4.34b) 

(4.34c) 


(Schwinger also explored the possibility of operator variations, see, for exam¬ 
ple, his les Houches lectures [Schwinger 1970a' .) We learn from the genera¬ 
tors, 

Gt = -H6t, Gg = J2PaSqa, (4.35) 


that the change in some function F of the dynamical variable is 


so we deduce 


SF = —St ■ 
dt 


-[F,Gl 


dF 

dt 

dF 

dqa 


\[F,Pa]. 


(4.36) 


(4.37a) 

(4.37b) 


Note that from this the canonical commutation relations follow, 

[qa,Pb] = iSab, [pa,Pb]=0, (4.38) 

as well as Newton’s law, 

1 f)M 

p, = —[H,pa] = -^. (4.39) 

* dqa 


If we had used Lp instead of Lq, we would have obtained the same equa¬ 
tions of motion, but in place of Gq, we would have obtained 


Gp — ^ ^ qgSpa, 


(4.40) 
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which implies 

dF 1 

= (4.41) 

From this can be deduced the remaining canonical commutator, 


[qa,qt]=0, (4.42) 

as well as the remaining Hamilton equation, 

1 f)M 

q, = -[q,,H] = ^. (4.43) 

It is easy to show that the effect of changing the Lagrangian by a total time 
derivative (which is what is done in passing from Lg to Lp) is to change the 
generators. 

We now turn to examples. 


4.1 Harmonic Oscillator 


The harmonic oscillator is defined in terms of creation and annihilation op- 
eratorsH and y, and the corresponding Hamiltonian il, 

[ 2 /,y'l'] = l, (4.44a) 

H = uj . (4.44b) 

The equations of motion are 

^ = -[y,H] = -ujy, (4.45a) 

at I i 

^ = = (4.45b) 

at I i 

Eigenstates of y and y^ exist, as right and left vectors, respectively, 

y\y') = y'\y'), (4.46a) 

{y%^ = y^'iy^'l (4.46b) 


^ We follow Schwinger’s usage of y for the annihilation operator, instead of the more 
usual a. 
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while (y'l and \y^') do not existH These are the famous “coherent states,” 
to whom the name Roy Glauber [Glauber 1963) is invariably attached, 
although they were discovered by Erwin Schrodinger jSchrddinger 1926] , 
and Glauber’s approach, as he acknowledged, followed that of his mentor, 
Schwinger [Schwinger 1953 . 

The transformation function we seek is therefore 

{y^\h\y",t 2 ). (4.48) 

If we regard y as a “coordinate,” the corresponding “momentum” is iy^: 

■ I dH . . dH 

y = -ujy=——^ lyf = -Ujyf = -(4.49) 

I my' oy 

The corresponding Lagrangian is therefore^ 

L = iy\y — H. (4-51) 

Because we use y as our state variable at the initial time, and at the final 

time, we must exploit our freedom to redefine our generators to write 

Wi 2 = j dtL - iy\ti).y{ti). (4-52) 

Then the variation of the action is 

SWi 2 = -iS(yl.yi) + Gi - G 2 

= -iSyl-yi - iyl-Syi + iyl-Syi - iy 2 - 5 y 2 - H Sti + H St 2 
= -iSyl-yi - iy\. 5 y 2 - H{Sti - 6 t 2 ). (4.53) 


Then the quantum action principle says 

= i{y^',ti\-iSyl'yi-iylSy2-u;yl'yi{Sti-6t2)\y",t2), ( 4 . 54 ) 
since by assumption the variations in the dynamical variables are numerical: 


® If {y'\y = y'{y'\ then we would have an evident contradiction: 


1 = {y'\ly,y'’]\y'} = y'iy'ly'^ly'} - {y'\y''^\y')y' = 0 . 


(4.47) 


^ We might note that in terms of (dimensionless) position and momentum operators 

iy^-y = - ip)-{y + w) = \{p-y - y-p) + t 4 :( 9 ^ +p^). (4-50) 

2 2 4 at 

where the first term in the final form is the average of the Legendre transforms in Lq 
and Lp. 
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[Syl,yi] = [yl,Sy2], 


(4.55) 


and we have dropped the zero-point energy. Now use the equations of motion 
(|4.45a|) and (I4.45bp to deduce that 




(4.56) 


and hence 


-iu;y^'e-^^^*^-*^\Sh-St2)y"\y",t2) 


= {y^',ti\y",t2)S 




(4.57) 


From this we can deduce that the transformation function has the exponential 
form 


{y^',ti\y",t 2 ) = exp 




(4.58) 


which has the correct boundary condition at ti = ^ 2 ! and in particular, (0|0) = 

1 . 

On the other hand, 


{y^',ti\y",t 2 ) = {y'^'le 


-iff(ii-t2)| // 


|y"), 


(4.59) 


where both states are expressed at the common time ^ 2 , so, upon inserting a 
complete set of energy eigenstates, we obtain (t = ti — ^ 2 ) 


Y,{y^'\E)e-^^^{E\y"), 


(4.60) 


which we compare to the Taylor expansion of the previous formula. 


n—0 


This gives all the eigenvectors and eigenvalues: 

En = nuj, n = 0,1, 2,..., 

{y^T 


{y^'\Er.) = 


{Eu\y") = 


n\ 

.ff\n 


{y”) 


(4.61) 

(4.62a) 

(4.62b) 

(4.62c) 


These correspond to the usual construction of the eigenstates from the ground 
state: 

K) = ^|0). 


(4.63) 
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4.2 Forced Harmonic Oscillator 


Now we add a driving term to the Hamiltonian, 

H = ujy^y+ yK*{t)+ y^K{t), (4.64) 

where K{t) is an external force {Kraft is force in German). The equation of 
motion is 

= wy + Nr(t), (4.65) 

while y^ satisfies the adjoint equation. In the presence of K{t), we wish to 
compute the transformation function {y"'^',ti\y'',t 2 )^■ 

Consider a variation of K. According to the action principle 

SK{y''',ti\y'',t2)^ = {y''',ti\i6KWi2\y",t2)^ 

=f dt[SKy'’+SK*y]\y",t 2 )^. (4.66) 
Jt2 

We can solve this differential equation by noting that the equation of motion 
(14.651) can be rewritten as 

4 ( 4 - 67 ) 


which is integrated to read 



e“*y(t) - e*‘"‘= 2 /(t 2 ) = A f dt' e^‘^*'K{t'), 

Jt2 

(4.68) 

or 

y(t) = i / dt'e-“(‘-*')A:(t'), 

Jt2 

(4.69) 


and the adjointd 


® The consistency of these two equations follows from 


AujtA , Aujto , 

e lyi = e ^y 2 — 


i dt’ K{t’), 
J to. 


(4.70) 


so that the adjoint of Eq. (14.6911 is 


[y(t)]t = e*-* 




i dt'i 
J to 


-icjt T^*(+l 


K*{t') 




+ i r dt' 

Jt, 


+ i r dt' 

J t2 

(4.71) 


which is Eq. (14.7211 . 
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/■ti 

yt(i) = -ij dt' (4.72) 

Thus our differential equation (14.661) reads 

MWM)^ J , / t/ ^ I ^ 

I:/,*.) 


= —i dt6K{t) 
Jt 2 

-i r dt5K*{t) 


-Mt-t2)y»_- f dt' 

Ju 


(4.73) 


Notice that in the terms bilinear in K and K *, K always occurs earlier than 
K*. Therefore, these terms can be combined to read 

fti , 

-Sk dtdt' (4.74) 
Jt 2 

where the step function is 



1, t > 0, 

0, t < 0. 


(4.75) 


Since we already know the K = 0 value from Eq. (14.581) . we may now imme¬ 
diately integrate our differential equation: 


{y^',h\y'',t 2 )^ = exp 


yP^-iu,(t^-t2)yn 


Jtn Ju 


t 2 

ftl 

- / dtdt' 

Ju 

The ground state is defined by y" = = 0, so 

nOO 

(0, fi|0, ^ 2 )^ = exp — dtdt'K*{t)r]{t — ^K{t') 

J —00 


(4.76) 


(4.77) 


where we now suppose that the forces turn off at the initial and final times, 
t 2 and ti , respectively. 

A check of this result is obtained by computing the probability of the 
system remaining in the ground state: 
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1(0, ti|0, ^ 2 )^ 1 ^ = exp I — y dtdt' K*{t)e 


= exp 
= exp [-\K{uj)\^] 
where the Fourier transform of the force is 


/ OO 

dtdt' 

-00 


K{uj)= / dte'^'K{t). 


(4.78) 


(4.79) 


The probability requirement 


|(0,tl|0,t2)^|^<l 


(4.80) 


is thus satisfied. We see here a resonance effect: If the oscillator is driven 
close to its natural frequency, so K{uj) is large, there is a large probability of 
finding the system in an excited state, and therefore of not remaining in the 
ground state. Let us calculate this transition amplitude to an excited state. 
By setting y” = 0 in Eq. (14.761) we obtain 


{y^',ti\0,t2)^ 


exp 


/ OO 

dte-^‘^(t,-t)K(^t) 

-OO 


(0,tl|0,t2)^ 


ti|n, ti){n, ti\0, fa)*”, (4.81) 


where we have inserted a sum over a complete set of energy eigenstates, which 
possess the amplitude [see Eq. (I4.62bp ] 


(yt» 



(4.82) 


If we expand the first line of Eq. (14.811) in powers of y^', we find 

I'—?')" 

(n,fi|0,t2)'^ = (4.83) 

Vn! 


The corresponding probability is 


p{n,0)^ = |(n,ti| 0 ,t 2 )*’|^ 




\K(u,)\^ 

5 


(4.84) 
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which is a Poisson distributioi:!^ with mean n = \K(uj)\'^. 

Finally, let us define the Green’s function for this problem by 

G{t - t') = 

It satisfies the differential equation 

as it must because [see Eq. (I4.65p ] 


- w ) y(t) = K(t), 


where y{t) is given by [see Eq. (I4.69p ] 

/ OO 

dt'G{t-t')K{t'). 

-OO 

Similarly, from Eq. (j4.72p 

/ OO 

dt' Git' 

-OO 

We can now write the ground-state persistence amplitude (14.12411 as 

noo 

(0,ti|0,t2)'^ = exp —i / dtdt'K*{t)G{t — t')K{t') 

J —OO 

and the general amplitude (14.761) as 


(4.86) 


(4.87) 


(4.88) 


(4.89) 


(4.90) 


(4.91) 


(y'f',ti|?/",t 2 )^ = exp I - i f dtdt'[K*{t) + iy^'6{t-ti)] 

L —OO 

xG{t - t') [K{t') + iy"6{t' - t 2 )] |, (4.92) 

which demonstrates that knowledge of (0,ti|0,t2)^ for all K determines ev¬ 
erything: 

{y^',ti\y",t 2 )^ = (4.93) 


® A Poisson probability distribution has the form pin) = A^e ^/n\. The mean value 
of n for this distribution is 


°° °° °° 

^ = E ’"pW = E )—WT = ^ E pW = 

n = 0 71 = 0 ^ n=0 


(4.85) 
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4.3 Feynman Path Integral Formulation 


Although much more familiar, the path integral formulation of quantum me¬ 
chanics [Feynman 1942[ [Feynman 1949[ [Feynman 1965[ is rather vaguely de¬ 
fined. We will here provide a formal “derivation” based on the Schwinger 
principle, in the harmonic oscillator context. 

Consider a forced oscillator, defined by the Lagrangian (note in this sec¬ 
tion, H does not include the source terms) 


L = iyKy - H{y,y^) - Ky'^ - K*y. 

As in the preceding section, the action principle says 

rti 

SK{ 0 ,ti\ 0 ,t 2 )^ = -i{0,ti\ I dt[SKy^ + SK*y]\ 0 ,t 2 )^, 
or for t 2 < t <ti, 


't-2 




(4.94) 


(4.95) 


(4.96a) 

(4.96b) 


where we have introduced the concept of the functional derivative. The equa¬ 
tion of motion 


iy - - K = 0 -iw'l' - — - K* = 0 

ly A u, ly A u, 

is thus equivalent to the functional differential equation. 


0 = <^ i 


K{t),W 




-iF(t) (0,ti|0,t2)^, 


where (the square brackets indicate functional dependence) 


W[y,y'’] = [ dt[iy\t).y{t) - H{y{t),y'<{t))]. 
Jt2 

The reason Eq. (14.9811 holds is that by definition 

(5 


6 K{t) 




(4.97) 


(4.98) 


(4.99) 


(4.100) 


so 
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*[ ^ ’4 \Jl^'M5K*{t') ~ ^\5K*{t'yJK{¥) 

. d i6 d fj ( *<5 \ 

^ ''^dK*{t) ~ d{iS/SK{t)) \dK*{t)' 5K{t)) ’ 


(4.101) 


which corresponds to the first two terms in the equation of motion (I4.97L 
under the correspondence 


y ^ ^ 


y ^^sK- 


(4.102) 


Since [[if, M^], Iff] = 0, we can write the functional equation (I4.98|) as 

0 = ^iW[iS/SK\^S/SK]J^^-iW[^S/SK\^S/SK]^Q^^^^Q^^^^^K_ (4.103) 


The above equation has a solution (up to a constant), because both equations 
(j4.97|) must hold, 

(0,ti|0,t2)^ = (4.104) 

where S[K], S[K*] are functional delta functions. The latter have functional 
Fourier decompositions (up to a multiplicative constant), 

6[K] = J (4.105a) 

S[K*] = J (4.105b) 

where [dy] represents an element of integration over all (numerical-valued) 
functions y(t), and so we finally have 

(0,ti|0,t2)^’^* 

]exp^-z^ dt[K{t)y^t) + K*{t)y{t)] +iW[y,y^] 

[dy][dy^]exp dt [iy^y - H{y,y^) - Ky^ - K*y]^ , (4.106) 

where y, y^ are now numerical, and the functional integration is over all 
possible functions, over all possible “paths.” Of course, the classical paths, 
the ones for which W — J dt{Ky"^ -\- K*y) is an extremum, receive the greatest 
weight, at least in the classical limit, where ?i —>■ 0. 
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4.3.1 Example 


Consider the harmonic oscillator Hamiltonian, H = ujy^y. Suppose we wish to 
calculate, once again, the ground state persistence amplitude, (0, ti|0, t 2 )^- 
It is perhaps easiest to perform a Fourier transform, 


y{t^) 



dt^^*y{t), 




dte-^^*y^{t). 


(4.107) 


Then 


/ °° Hu 

dtyHt)y{t) = J ^y{v)y*{-u), 

/ OO poo J 

dtiy\t)y{t) = j —uy{u)y*{-u). 
-00 J — OO 


(4.108a) 

(4.108b) 


Thus Eq. (14.1061) becomes 

(0,ti|0,t2)^’^* 

= J[dy][dy*]exp^^i J ^[y{u){u - uj)y*{-u) 


- y*{-u)K{u) - y{u)K*{-u)] 

K{v) 


yiv) - 


= J [dy][d 2 /*]exp|i J ^ 

J ZTT U — LU 

= J[dy][dy*]expi^i J ^y{v){u - uj)y*{-u) 


{v-uj) 


y*{-i^) - 


K*{-u) 


(4.109) 


since the functional integral in the third equality, obtained by shifting the 
integration variable, 


y{v) - y{v), 


y*{-v) - 


U — LO 

K*{-u) 


y*i-^), 


(4.110a) 

(4.110b) 


is (0, ti|0, t 2 )^“'^ = 1- How do we interpret the singularity at = w in 

the remaining integral? We should have inserted a convergence factor in the 
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original functional integral: 

exp J ^ [... ]^ ^ exp J ^ [• • • + iey{v)y*{-iy)]^ , (4.111) 

where e goes to zero through positive values. Thus we have, in effect, v — oj^ 
V — uj + ie and so we have for the ground-state persistence amplitude 

(0,<i|0,t2)^’^* = (4.112) 


which has the form of Eq. (14.911) , with 



dv 

27r iz — w -|- te ’ 


(4.113) 


which is evaluated by closing the v contour in the upper half plane if t—t' < 0, 
and in the lower half plane when t — t' > 0. Since the pole is in the lower half 
plane we get 

G{t - t') = -ir]{t - (4.114) 

which is exactly what we found in Eq. (j4.86|) . 

Now, let us rewrite the path integral ()4.10611 in terms of coordinates and 
momenta: 


q = 




y = 




(4.115a) 

(4.115b) 


Then the numerical Lagrangian appearing in (14.1061) may be rewritten as (see 
footnote m above) 


L = iy^y — coy^y — Ky^ — K*y 


■ W / -P\ l ■ , ■P\ I 2 , P 


UJ 






— '^Kq — —^Kp 

(jJ 


= -^w + L{q,q,t), 


(4.116) 


where, if we set q= p, the Lagrangian is 
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L{q,q,t) 



1 


+ Fq, 


if 


= 0, F = -V^^K. 


In the path integral 


[dy][dy^] = [dq][dp] 


d{q,p) 


where the Jacobian is 


5(2/, y^) 




d{q,p) 


i 

\/2u) 

i 

\/ 2a; 


(4.117) 

(4.118) 


(4.119) 


(4.120) 


and so from the penultimate line of Eq. (|4.116p . the path integral (j4.106|) 
becomes 


(0,ti|0,t2)^ 


J [dy][dy^ exp 


dtL{y,y^ 


't2 


j [dq][dp]exp 





(4.121) 


Now we can carry out the p integration, since it is Gaussian: 

j [dp]Fd'dt[-hp^+pq] ^ J ^dp]F f dt[-i{p-qf + ^q^] 

/ oo 

(4.122) 

-OO 


Here we have discretized time so that p{ti) = pi, so the final functional 
integral over p is just an infinite product of constants, each one of which 
equals At. Thus we arrive at the form originally written down by 

Feynman [Feynman 1965| , 

(0,ti|0,t2)^ = y[dg]exp dtL(g,g,t)|, (4.123) 


with the Lagrangian given by Eq. (I4.117p . where an infinite normalization 
constant has been absorbed into the measure. 
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4.4 Toward Source Theory 

Let us return to the action principle. Recall from Eq. (14.1121) 

(OtilOtz)^ = (4.124) 

The action principle says 

S{h\t2) = i{h\S[Wi = J dtL]\t2). (4.125) 

In a general sense, the exponent in Eq. (I4.124p is an integrated form of the 
action. In solving the equation of motion, we found in Eq. (14.891) 

y(t) = + J dt'G{t - (4.126) 

where the first term is effectively zero here. The net effect is to replace an 
operator by a number: 


y'{t) = J dt'G{t-t')K{t'). (4.127) 

Then Eq. (14.1241) can be written as 

(OtilOta)^ = e-* W. (4.128) 

Recall that the action was was the integral of the Lagrangian (14.941) . or 


W = dt 


- wy'l'y - y^Kit) - yK*{t) 


(4.129) 


so we see one term in Eq. (14.1281) here, and the equation of motion (14.651) 
cancels out the rest! So let’s add something which gives the equation for y': 

(OtilOts)^ = e^fdt[y^'i4!:y'-‘^y*’y'-y^'K-y'K-‘] ^ ^^w^ (4.130) 

Now insist that W is stationary with respect to variations of y', y^', and we 
recover the equation of motion, 

y'W =(4.131) 


This is the starting point for the development of source theory, which will be 
treated in Chap. [7l 



















Chapter 5 

Time-cycle or Schwinger-Keldysh 
formulation 


A further utility of the action principle is the time-cycle or Schwinger- 
Keldysh formalism, which allows one to calculate matrix elements and con¬ 
sider nonequilibrium systems. Schwinger’s original work on this was his fa¬ 
mous paper [Schwinger 1961| ; Keldysh’s paper appeared three years later 
[Keldysh 1964| , and, rather mysteriously, cites the Martin-Schwinger equilib¬ 
rium paper [Martin 1959] . but not the nonequilibrium one [Schwinger 196l[ . 
The following was extracted from notes from Schwinger’s lectures given in 
1968 at Harvard, as taken by the author. 

Consider the expectation value of some physical property F(t) at a par¬ 
ticular time ti in a state [ 6 ,^ 2 ): 

= ^( 6 't 2 |a'ti)(a'|F|a")(a"ti| 6 't 2 ), (5.1) 

a'a" 

which expresses the expectation value in terms of the matrix elements of the 
operator F in a complete set of states defined at time F, {|a'ti)}. Suppose 
the operator F has no explicit time dependence. Then we can use the action 
principle to write 


S{a'ti\b't2) 


i{a'ti\6 



\b%), 


(5.2a) 


and so 


S{b't2\ati) 


—i(b't2\S 



(5.2b) 


which can be obtained from the first equation by merely exchanging labels, 



(5.3) 


If we consider 
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{b't2\b't2) = '^{b't2\a'ti){a'ti\b't2), (5.4) 

a' 

the above variational equations indeed asserts that 

S{b't2\b't2) = 0. (5.5) 

We can interpret the above as a cycle in time, going from time ^2 to ti and 
then back again, as shown in Fig. 15.11 But, now imagine that the dynamics is 



h 


Fig. 5.1 A “time-cycle,” in which a system advances forward in time from time t 2 
to time ti under the influence of a Lagrangian L+, and then backward in time from 
time ti back to time t 2 under the influence of Lagrangian L_. 


different on the forward and return trips, described by different Lagrangians 
L+ and L-. Then 


S{b't2\b't2) 


i{b't2\S 


r 

/•ii 1 

/ dtLj^ — 

/ dtL_ 

Jt2 

It2 


\b't2). 


In particular, consider a perturbation of the form, 


(5.6) 


H = Ho + X{t)F, (5.7) 

where X(t) is some time-varying parameter. If we have an infinitesimal change, 
and, for example, (5A+ ^ 0, SX- = 0, then 

rti 

Sx^{b't2\b't2)^+^- =-t{b't2\ dtSX+F\b't2). (5.8) 

J t2 


If we choose i5A+ to be an impulse, 
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in this way we obtain the expectation value of F{t'). 

Let’s illustrate this with a driven harmonic oscillator, as described by 
Eq. (14.641) . so now 


H+= + Kl{t)y + K+{t)y\ (5.10a) 

H- = Loy^y + K^{t)y-\-K-{t)y\ (5.10b) 

which describes the oscillator evolving forward in time from t 2 to ti under 
the influence of the force and backward in time from ti to t 2 under the 
influence of AT-, as shown in Fig. 15.21 From the variational principle we can 



Fig. 5.2 A time cycle In which the harmonic oscillator evolves from time t 2 to time 
ti under the influence of a force K^, and then from ti back to time t 2 under a force 


A.. 


learn all about y and yK We have already solved this problem by a more 
laborious method above, in Section 

It suffices to solve this problem with initial and final ground states. If we 
consider only a AT* variation, 

rti 

SK-{0t2\0t2)^+’^- = -i{0t2\ / dt [6Kl{t)y+{t) - SK*_{t)y_{t)] |0t2). 

Jt2 

(5.11) 

Now we must solve the equations of motion, so since effectively 1 /(^ 2 ) 0, 

we have from Eq. (14.691) . 


y+it) = -i f dt' 
J t2 

pti 

y-it) = -i dt' 

J t2 


(5.12a) 



dt'e-*“(‘-‘')A:_(f').(5.12b) 


The last term in the second equation is 
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fti 

Jt2 


(5.13) 


so naming the advanced and retarded Green’s functions by extending the 
definition in Eq. (14.861) . 


G, 


I , (5.14) 


which satisfy the same differential equation (I4.87L we effectively have 

't2 

y 


y+{t) = j dt' Grit — t')K+(t'), (5.15a) 

Jt2 

_it) = -i f\t'e-^‘^^*-*'^K+it)+ [" dt'Gait-t')K_it'),i5.15h) 
Jtn Jto 


t2 ^t2 

The solution to the variational equation f|5.1ip is now 

^0^210^2) 

(5.16) 


This should reduce to 1 when K+ = K- = K, so 

- iGrit - t') + iGait - t') + = 0, (5.17) 

which is, indeed, true. 

As an example, consider K-it) = Kit), K+it) = Kit + T), that is, the 
second source is displaced forward by a time T. This is sketched in Fig. 15.31 
What does this mean? From a causal analysis, in terms of energy eigenstates, 
reading from right to left, 

(0t2|0t2)*’-’^+ = (5.18) 

n 

The effect on the second transformation function is the same as moving the 
n, tl state to a later time, 

(nti|0t2)^^‘+^^ = (nil + r|0t2)^^*^ = e-“‘"^(nti |0t2)^^*^ (5.19) 


so this says that 


(0t2|0t2)'^-'^+ =^e-““^p(n,0)^. 


n 


(5.20) 
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Fig. 5.3 Time cycle in which K-{t) = K{t), K+{t) = K{t + T), that is, the forces 
are the same on the two legs, but displaced in time. 


which gives us the probabilities directly. From the formula (15.161) we have, 
using Eq. (I5.17p . 

= (5.21) 

where 

j = J dte^‘^*^K{t). (5.22) 

Thus we immediately obtain Eq. (14.841) . or 

p(n,0)^ = (5-23) 

nl 

The above Eq. (15.211) can be directly used to find certain average values. 
For example, 

(e—= (5.24) 

Expand this for small ojT and we find 

(n)^ = ItP- (5.25) 

In a bit more systematic way we obtain the dispersion: 

_ g|7|^(e“’“'^-l-|-iwT)_ (5.26) 


Expanding this to second order in ojT we get 















56 


5 Time-cycle or Schwinger-Keldysh formulation 


((n - (n))2) = (n^) - (n)^ ^ {Anf = I 7 P = (n), (5.27) 


or 


An 1 
(n) 


(5.28) 


For large quantum numbers, which corresponds to the classical limit, the 
fluctuations become relatively small. 

Now consider a more general variational statement than in Eq. (I5.11|) , 


5( I )^-^+ =-*(■/ dt[6Kly+ - 6K*_y. + 6K+yl - 6K.yl\ , (5.29) 


we see that since we can change the source functions at will, and make very 
localized changes, it makes sense to define the variational derivatives 

I )^* = ( \y+(.t)\ (5.30a) 

= (5.30b) 

All expectation values of operator products at any time can be obtained in 
this way—in particular, correlation functions. Repeating this operation we 
get 


(-*) 




(<2|»+(i')|<2) 


K± 


= {t 2 \yli.t)y+ ■ (5.31) 


The operators are multiplied in the order of the time development. The only 
place where K- appears is in the latter part of the time development. See 

Fig. [El 

The distinction between ± disappears if we now set A"+ = Ar_: 


SK_{t) 




= {0h\yHt)yit')\0t2)^. (5.32) 


K+=K-=K 


As an example, set t = t' = ti] then, from Eq. (I5.16|) . this reads for the 
number operator N{t) = y^t)y{t), 




j dtK*{t)Ga{t-ti) 


dt'Griti - t')K{t') 





|7P, 

(5.33) 


as before, Eq. (15.251) . 
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Fig. 5.4 Variational derivatives pick out operators at definite times t and t'. 


We would like to use more general starting and ending states than the 
ground state. We can obtain these by use of impulsive forces. It is convenient 
to deal with all states at once, as in the generating function for p(n, 0 )^ 
considered above. Think of a time cycle starting at time t 2 , advancing forward 
to time ti, during which time the force acts, then moving back in time to 
a time under the influence of the force K -—See Fig. 15.51 Let ^2 = ^ 2 + t. 



Fig. 5.5 Time cycle with different forces, and K-. on the forward and backward 
moving segments. Now the initial time of the time cycle, t 2 , is different from the final 
time of the time cycle, tj, with r = — t 2 . It is assumed that the time ti is later 

than both t 2 and tjj and that the forces are localized as shown. 


This displacement injects energy information. Consider 

Y,{nt2\nt2)^^ = tr(t'2|t2)^± = ^ , (5.34) 


n 


n 
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which uses (no force acts between times ^ 2 ) 




(5.35) 


Analysis of this formula will yield individual transformation functions. 

Now we must solve the dynamical equations subject to boundary condi¬ 
tions. Let us compare tr(t 2 |j/-i-(t 2 )|t 2 ) with tr(t 2 |y-(t 2 )l^ 2 )- The first is 

ta'(^ 2 l 2 /+fe)|i 2 ) = '^{nt 2 \y+{t 2 )\nt 2 ) = '^{nt 2 \nt 2 ){n\y\n), (5.36a) 

n nn' 

while the second appears as 

ti'(^2ly-(4)l^2) = ^{n't2\y-it2)\n't2) = ^{n'\y\n){nt2\n't2). (5.36b) 

n' nn' 


Here, by introducing a complete set of states at the time of the operator, 
we have expressed the formula in terms of the matrix elements of stationary 
operators. Remarkably, we see that the two expressions are equal; in effect, 
there is a periodicity present here: 


y+ih) =y-it2), (5.37) 

as far as traces are concerned. Now, the equations of motion (|4.65ll for the 
operators read 

- w) y(t) = Ar(t), (5.38) 

which has solution (|5.12bl) with the addition of the initial term, or 


y-{t) = e-“(*-‘=)y+(t2) -i dt' A:+(t') 

J ±2 

pti ^ 

+ i J dt'e-*“(*-‘)A:_(t')- (5.39) 

In particular, 

y.it’2) = e-^-^y+{t 2 ) - IJ - K_){t'). (5.40) 


Note that the integrals sweep over the full force history. Let us let t 2 = 0 
for simplicity, although we will keep the label. Because of the periodicity 
condition (15.371) this reads 


or 


- 1) y+{t2) = -* y dte^‘^\K+ - K_){t) = - 1 ( 7 + - 7 _), (5.41) 
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y+(^ 2 ) = - 7-)- 

What we are interested in is 

tr(^ 2 l^ 2 )^" 

tr(t^|t 2 ) 


(5.42) 


(5.43) 


The denominator, which refers to the free harmonic oscillator, is immediately 
evaluated as 

OO 

= = (5.44) 

n—0 

(If T be imaginary, we have thermodynamic utility.) We have then the vari¬ 
ational equation 




’tr(tyt 2 ^’ 

_ tr(t' 2 |t 2 ) . 


I dt {5Kly+ - 5K*_y_) [^ 2 )^=^ 
tr(t' 2 |t 2 ) 


(5.45) 


Exactly as before, we get an equation for the logarithm—looking at the 
previous calculation leading to Eq. (I5.16|) , we see an additional term, referring 
to the y+{i 2 ) boundary term in Eq. (I5.40|) . The periodic boundary condition 
then gives 

- ^ ^^(7+ - 7-)(7-f - 7-)- (5-46) 

Therefore, to convert (0t2|0t2)^* in Eq. (15.161) to 


ti'(t2|i2) ^ ■ ’ 

we must multiply by 


This holds identically in r; in particular, in the limit where r —> —loo, which 
corresponds to absolute zero temperature, we recover ( 0 ^ 210 ^ 2 )^"*" ■ 

We find, generalizing Eq. (I5.16p 

^ ^-^fdtdt'Kl{t)GAt-t')K+(t') 

g —mCJT 

^ gZ Jdtdt'K'^{t)Go{t — t')K^{t')^J dt dt' 


which is the exponential of a bilinear structure. This is a generating function 
for the amplitudes {nt 2 \nt 2 )^^ ■ But it is useful as it stands. 
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Put r = —*/3; then this describes a thermodynamic average over a thermal 
mixture at temperature T, where /3 = 1/kT in terms of Boltzmann’s constant, 





In terms of this replacement, 


1 

^iuJT _ 1 



(5.50) 


(5.51) 


because 

- = - - -ln(y^ = - - -In- ^ - = - - - 

^^g-zna;r d{-iUJT) ' d{-iujT) 1 - - 1 ' 

(5.52) 

Now consider a time cycle with displacement T: the system evolves from 
time t 2 to time ti under the influence of the force K^(t), and backwards in 
time from ti to ^2 under the force K-{t): 


K_{t)=K{t), K+{t) = K{t + T). (5.53) 

This is again as illustrated in Fig. 15.51 with these replacements. What is the 
physical meaning of this? Insert in Eq. (15.491) a complete set of states at time 
ti- 

{nt2\nt2)^^ = ^(nt2|n'ti)'^~(n'ti|nt2)^+- (5.54) 

n' 

We did this before for the ground state. The effect is the same as moving the 
starting and ending times. Appearing here is 

{n'h\nt 2 )^^*+^^ = (n'ti + r|nt 2 + = e-“'‘"^(n'ti|nt 2 )^(‘)e“‘"^. 

(5.55) 

Therefore, 

n)^ = (5.56) 


Therefore, as a generalization for finite r of Eq. (j5.26|) . we have from 
Eq. ()5.49p 


-1 


E' 


^ ^ ^ — inujT — —n)ujT'^l 


= exp 


(e-“" - 1) |7P - (e-" - 1) (e-- - l) I 7 I 


,(5.57) 


where T gives the final state, and r the initial state. This used the observation 
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/' 


Expand both sides of Eq. (15.571) in powers of T, and we learn 

p—inujT 

- tujT Y,{N - n)^ = -iwTlqp, 

Z^n' 


(5.58) 


(5.59) 


or 

(iV-n)f = |7|^ (5.60) 

which generalizes the earlier result (j5.25l) . Now apply Eq. (15.591) as a gener¬ 
ating function, 

(iV-n)^ = |7P, (5.61) 

which reflects the linear nature of the system. 

We can rewrite the above generating function more conveniently, by mul¬ 
tiplying by 

(5.62) 

that is, Eq, (15.571) can be written as 


1 


Ee- 


= exp 


^ ^ ^—inujT —{N))ijjT\^K 


(e-*‘^^-l + iwT)|7|2- 


g*a;r _ I 


(e-“^ - 1 ) - 1 ) | 7 | 


(5.63) 


Now pick off the coefficient of —(u;r)^/2: 

(«))»)„* = + (5.64) 

or 

((iV-(iV))2)^ = |7ni + 2(n)^]. (5.65) 

If, instead, we multiply Eq. (15.641) through by ^ ^gg ^j^jg gg 

a generating function, and learn from Eq. (I5.52|l that 

((iV-(iV))2)^ = |7p(l + 2n). (5.66) 

Note the simplicity of the derivation of this result, which does not involve 
complicated functions like Laguerre polynomials. 






















Chapter 6 

Relativistic Theory of Fields 


This section is an adaptation of Chapter V of lectures given at Stanford by 
Julian Schwinger in 1956 [Schwinger 1956| . 

A state of a physical system is defined in terms of the maximum number of 
compatible measurements which can be made upon the system. If the state 
were defined on a space-like surface (one in which all points are in space-like 
relation: (Ax)'^ — {At)‘^ > 0) then a measurement at any point is compatible 
with one made at any other point, since the disturbances introduced by the 
measurements cannot propagate faster than c, and hence cannot interfere. 
Thus, a state can be specified as an eigenvector of a complete set of com¬ 
muting, Hermitian operators a, associated with a definite space-like surface 
(T : \a',a). There always exists a coordinate system in which the space-like 
surface a is all of three-dimensional space at a given time; in this particu¬ 
lar Lorentz frame the state vector is just: |a',t). The problem of relativistic 
dynamics is to find the transformation function 


1^2 ^2) ■ 


( 6 . 1 ) 


a + da 

I I 

a 


Fig. 6.1 Spacetime volume bounded by two spacelike surfaces, cr and cr-t-dcr. Points 
in the blue and green shaded regions cannot interfere. 
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6 Relativistic Theory of Fields 


As in the non-relativistic case, we assume the existence of an action oper¬ 
ator Wi 2 such that 


S(a[(7ila2a2} = i{ai(Ji\5\Wi2]\a2<y2) ■ ( 6 . 2 ) 

The contributions to the action operator are now given by 

W^2=Y.^-+d-,-- (6-3) 

(7 


Since measurements made at points in the space-like shell da cannot interfere 
(e.g., in the two shaded areas, shown in Fig. 16.1|) their contribution to W 12 
is additive, 

W 12 = f {dx)C{x) = f {dx)C[xa{x)], (6.4) 

J (72 <72 

where the Xa{x) = Xa{x,y, z,t) are the dynamical variables of the system, 
necessarily Hermitian operators; the relativistic requirements automatically 
introduce the concept of fields. The relativistic notation used has the form 


x^ = t = —Xq, x^ = Xfc, where fc = 1, 2, 3, 

d 

{dx) = dx^dxidx 2 dx 3 , 
where the metric used is 


/-I 0 0 0\ 
0 10 0 
0 0 10 
\ 0 0 0 1 / 


(6.5) 


( 6 . 6 ) 


The action principle again states that for a given dynamical system, the 
variations arise only from the end-point, that is, 


SW12 — Gi(tTi) — 02(^2) 


(6.7) 


As before, from this requirement follow the equations of motion and the 
generators of infinitesimal transformations which yield the commutation re¬ 
lations of the field operators. The Lagrangian density £ which will yield first 
order field equations is 

^ = ^ - df^x^^x) -'^ix), ( 6 . 8 ) 

where the are a set of four numerical matrices, and space and time deriva¬ 
tives appear on a symmetric basis. The symmetrization of the kinematical 
term relates to the possibility of adding to £ the relativistic analogue of our 
previous total time derivative, a four dimensional divergence. If 
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C = C- (6.9) 

then 

Wi2 = Wi2- f {dx)df,f'" =Wi2 - f + f ( 6 . 10 ) 

J (72 'J (y\ 'J (72 

and 

G(ai) = G(ai) - / da^r, G(a2) = G(a2) - [ da^f^. (6.11) 

*/ CTl *-> <72 

As before, the equation of motion is unchanged. Since £ is to be Hermitian, 
we require 

so n'f = n (6.12) 

The rank of that of the number of independent fields. Note that the 
variation 5{a'iai\a2cr2) is independent of any coordinate system, since £ is a 
Lorentz scalar. 

We can now infer some fundamental properties from the requirement of 
invariance of C. Consider the coordinate transformation (Poincare transfor¬ 
mation) 

- r (6.13) 

where 

= 5A«, (6.14) 

We can divide the transformations into two subsets by considering the effect 
of (16.131) on goo, 

— ffoo = 1 = (Co) ~ ^ (Co) (Co) = 1 + ^ (Co) >1- (6.15) 

k—1 k 

Since £°o = it follows that we can never make a continuous change from a 
positive to a negative sense of time, i.e., generate an improper transformation 
continuously. We shall consider only the group of continuous proper Lorentz 
transformations. 

Under such a coordinate change, the fields Xaix) change to new ones Xa^) 
connected by a real linear transformation: 

x{x)=Lx{x). (6.16) 

Using dfi = and writing the scalar C in the new frame we find 

\ {x^^d^X - d^^xA^x) - ^ (x) 

= i [x (LTyi^r^L) d,x - d.x (r^L^A'^L) x]-n (Lx) • (6.17) 
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Thus 

L^A^L = n {Lx) =n{x). (6.18) 

If we choose for Ti. the special form 

n = xBx = xL^BLx, (6.19) 

where B is Hermitian and non-singular, then = BL~^B~^. From the hrst 
of equations (16.181) . we then obtain 

L-^ {B-^A^^)L = ^^^{B-^A'') , ( 6 . 20 ) 

showing that the combination B~^A^^ transforms like a vector; this is implied 
by the required invariance of the kinematical term of C. 

Consider now a general inhomogeneous infinitesimal Lorentz transforma¬ 
tion, of the form 

= + ( 6 . 21 ) 

where the and are infinitesimals, and the rotational nature of the 
is expressed by the relations Then L can be written as 

L = 1 -I- where only the rotational terms appear, since the 

translations do not effect the gradient operators d^, and no corresponding 
changes in the Xol required to keep C invariant. The are (imaginary) 
operators, acting on the field variable, which will express the spin character 
of the helds. 

The variation SWa allows us to change the field components at each space- 
time point (call these changes SaXa), and to change the region of integration 
by displacing the boundary surfaces ci and tT2. In the previous non-relativistic 
treatment, instead of varying the end-point times and ^2 ,we used ^(t). 
Proceeding similarly here, we can express the variation of a space-like sur¬ 
face by varying the space-time coordinates under an infinitesimal Lorentz 
transformation (dcc^ = C/i ~ so chosen that on ui and (72 the required 

displacement is obtained. The change in the action is 
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SWi2 = i5o 


(dx) 


- (xA'^df,x - d^xA^x) - 


+ / dcTnSx^C — / danSox^C 


{dx) [SoxA^d^x - dfiXA^SoX - ^o'H] 




+ 


{dx) df, {xA>"Sox - SoxA>^x) 


dafiSx^C — 


da^dx^C 


£ {dx) 5^C + ^ da^ i (x^'^^oX - <5oX^''x)| 

-/ rfcr^|'^a:'"^ + ^(x^^^oX-<^oX^^x)|- 


( 6 . 22 ) 


Applying the action principle, the interior variation 5qC must vanish, giv¬ 
ing the field equations of motion. What remains is the difference of two gen¬ 
erators, G(cri) — G(o- 2), where 


G{a) = ^ dcj^ f^Sx^C + i (xA'^JoX - ^oX^^x)} • (6. 


23) 


To re-write (16.231) we recognize that the total change in the fields is due to 
both the ^0 variation of the fields at a given space-time point on a, and 
also to the variation induced by the infinitesimal Lorentz transformation of 
coordinates as a is displaced to cr -I- da; the latter is obtained from 


X {x) = Lx{x) = x{x) + -e^uS^'^x{x), 


But 

and therefore 


X (a;) = x(a:) - Sxf^df,x{x), 


X{x) - x{x) = Sx^dfj,x{x) + -efj,„S^''x{x)- 


(6.24) 

(6.25) 

(6.26) 


The right hand side of (16.2611 is then the variation induced in the Xa by the 
coordinate transformation; the total variation of the fields is then 


Sx{x) = Sox{x) + Sx^d^x{x) + -e^^S^''x{x)- 


(6.27) 


Solving for Jqx, and substituting into (16.2311 . we obtain for the generator 
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G(cr) = J da^ i ixA'^Sx - SxA>^x)^ 

- i I {xA^d-'x - d^'xA^x) Sx, 

- {xA^^S^-'x - S^'^xA^x) | • (6.28) 


Adding a surface term, it is possible to bring G{a) into the form: 

G(a)= J da^^ixA>^6x-6xA>^x) + J da^dx,T>^\ 

where is the symmetric stress-tensor operator 







xAi‘'S^^^X- S^^^A^x 


(6.29) 


(6.30) 


and the brackets { } represent symmetrization with respect to /t and ly. 

Applying the stationary action principle to the variation Sxi,, taken as 
arbitrary, we note that the invariance of the action operator implies the con¬ 
servation law: 


da^T'^'' = 


da^T^’', 


(6.31) 


which, in turn, implies the corresponding differential conservation law 


= 0. (6.32) 

The generator (16.291) can be split into two parts, one representing changes 

induced by the coordinate variation, and the other giving the variation in¬ 
duced by a change in the field variables, 

Gx= J daj^ {xA^dx - dxA^^x) , (6.33a) 

Gx = j da^T^’'5x, = (6.33b) 

where 

P'' = J da^T^'^, J^'' = j dax [x^^T^'' - ; (6.34) 

P" and are the generators for translations and rotations, respectively, and 
their commutation relations are determined by the group of transformations 
they represent. Specifically, P'^ is recognized as the 4-momentum operator, 
and as the relativistic generalization of the angular momentum operator. 
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The field equations are obtained by the vanishing of in (I6.22L 

2A>^d^,x = ^ or - 2A>^^d^x = (6.35) 

in terms of left and right derivatives, corresponding to the two equivalent 
ways of writing G^: 

Gx = j da^,xA'"^X, or = - j da^SxA^^x- (6-36) 

If we continued with these two pairs of expressions,we would obtain two 
forms for the commutation rules of y; their equivalence then leads to the 
requirement that the and x must decompose: 

X = 4> + 'il’, A^=a^ + s^, (6.37) 

where the are anti-symmetric and real, the are symmetric and imagi¬ 
nary, and the (j) and ip represent the kinematically independent fields of the 
Bose-Einstein and Fermi-Dirac types, respectively. The field equations are 
then 

Similarly, G^ can be. divided into its (p and ip parts, G^ = G^ + G^, 

G^ = j da^(pa^6(p = J da^ {a^5(p) cp, 

G^ = J da^ips^dip = J da^ {s^Sip) ip. (6.39) 

Again, we find (where { , } denotes an anticommutator) 

{Sipa{x),ip/ 3 {x'} = 0, 

[S(pa{x),(pk{x')] = [Scpk{x), 1 paix')] = [Slpa{x), (pk{x')] =0. (6.40) 

The first of equations (j6.40l) combined with the second of equations (j6.38bp 
implies that 'H must be an even function of ip. 

The field equations may be written as equations of motion by singling out 
the time differentiation, 


2A°dox 


dm 

dx 


2A'^dkX- 


(6.41) 
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If we took to be non-singular, we would be able to solve (16.411) for dOx- 
More generally, we now recognize the existence of the following possible sit¬ 
uations: 

1. is non-singular. In this case, all of the variables are kinematically inde¬ 
pendent. An example of this situation is the Dirac-Majorana spin-i field. 

2. AO is singular, but there are enough relations among the variables to de¬ 
termine all of them. Here, only those variables which possess equations of 
motion are kinematically independent, but the non-independent fields are 
determined from the independent fields. Examples of this are the spin zero 
and spin one fields. 

3. AO is singular, and there are not enough relations among the variables to 
determine all the fields, as in case [5] The classic example of this is the 
spin 1, zero mass, electromagnetic field, where the lack of determination 
corresponds to the possibility of introducing a gauge transformation. 


6.1 Inference of Particle Properties 


We now consider the generators of infinitesimal coordinate (Lorentz) trans¬ 
formations Gfc and the commutation relations they imply. From (I6.33bl) we 
have 

G, + (6.42) 

which, when applied to the space-time coordinates a;^, generates the new 

+ (6.43) 

which is Eq. (16.211) . Accompanying this transformation, we have the apparent 
change in the fields given by (16.261) . 


- Sxix) = x{x) - x(x) = Sx>^df,x(x) + ]^t^yS^''x{x), (6.44a) 

or 


- 5x{x) 








X{x), 


(6.44b) 


using = —Cvn- Comparing with [y, G^] = —iSx, for arbitrary translations 
(e^) and rotations (e^jj/), we obtain 


[X,P''] =-W'x: (6.45a) 

[y, = -i - x''d>^ + X- (6.45b) 


From (I6.45bl) the identification of with the intrinsic spin characteristics 
of the particle is evident. Considering the time component (I6.45al) we obtain 
the standard commutator equation of motion 
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-id°x = idox 



(6.46) 


To determine the manner in which the particle interpretation enters, we now 
represent the fields by Fourier integrals, 


Xi.x) = J (dp)e“'‘P^‘x(p), 


(6.47) 


where the xip) s-re operator functions of the numbers Substituting into 
(j6.45all and equating coefficients, we obtain 

[x{p),P''] = P^xip), (6.48a) 

or 

P^'xip) = xip) {P"' - p'") ■ (6.48b) 

Since the 4 operators P'^ all commute, we can have simultaneous eigenstates 
\P'), where P''\p') = P'"\p'). Applying (I6.48al) to these states 

P'^xW) = X iP'^' - pn \P') = iP'" - pI x\P'), (6.49) 

or writing x\P') some new eigenvector \p' — p), 

P'^IP' -p) = iP'" -pn\P'-p), (6.50) 

which shows that the effect of x on \P') is to produce a state whose momen¬ 
tum eigenvalue has been changed by —p'^, indicating the capacity of the field 
to absorb or emit p‘^. Specifically, considering the time component ^ = 0, we 
have two possible situations, depending on whether p° is greater than or less 
than zero. If p° > 0 , application of x on \P') yields a state with lower energy 
(pO' — p°'), and conversely if p° < 0. Since x is Hermitian, x can be split, in 
a Lorentz covariant way, into two parts, x = + X^~\ where the super¬ 

scripts indicate the positive and negative nature of the p^ terms which enter 
in the respective Fourier transforms, and where (x^"*"^)^ = X^~^. Then 
applied to a state of definite energy acts as an energy annihilation operator, 
and x^~^ 3.8 an energy creation operator. 

The x*'"''^ and x^~^ can then be used to create any physical state from the 
vacuum state, where we take the latter as the unique lowest-energy state. 
This vacuum state must necessarily correspond to the eigenvalue = 0, 
(fc = 1,2,3), since if one of the P^' were not zero, a rotation of the coordi¬ 
nate system could yield three non-zero momentum components, requiring (a 
super-position of) the corresponding eigenvectors for its description. But this 
precludes a description of the vacuum by a single non-degenerate state; we 
must therefore require that each P^' = 0. Since the vacuum is to have the 
lowest energy possible, and we may arbitrarily take this to be zero, we then 
characterize the vacuum state as that unique state for which P^^' = 0, where 
x(+) I 0 >= 0. 
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So far we have considered x as representing general fields; to introduce par¬ 
ticle properties consider the Fourier transform of x, and imagine the numbers 
Pfi related by the relation: 

— Pi^p^=m^, = ±\/ m? + p^. (6.51) 

Then for example, remains an operator which annihilates energy, but 

now is correlated with a momentum decrease. This is just the usual particle 
interpretation; if has this character it may be spoken of as a particle 
annihilation operator, and conversely for x*'”^- 


6.2 The Connection Between Spin and Statistics 


We consider the simplest system, corresponding to linear field equations. 
Taking % = X-^Xi this becomes for our two types of fields % = + 

where (j), ip represent fields of the I'** and 2”'^ kinds, respectively, 
meaning bosonic and fermionic fields. (Note that % must be even in ip, and 
therefore, as noted below Eq. (I6.40L no term of the form (pip can occur.) The 
matrix B (and therefore B^^'> and B^‘^'>) is Hermitian. Writing 

5^n = 6(PB''^'>(P + (PB^^'^dcp = 5(p (^B^^U + , (6.52) 


and since 


dm _ dm 

d(p d(p 


(6.53) 


for these variations of fields of the first kind, both terms in the above bracket 
must be equivalent, implying that B^^'> = Since B^^'> is also Hermitian, 

it is real. For Sip variations, the identical argument shows that B^^^ = 
i.e., is imaginary. It is precisely here that the connection between spin 
and statistics arises: We can construct matrices for B~^a^ only for particles 
of integral spin, and for B~^s^ only when the particle spin is an integer plus 
one-half. [See Eq. (j6.37l) .] 

Writing the equations of motion (16.3511 for either held 




(6.54) 


and assuming—as the simplest case—that B is non-singular, we obtain 


B ^A^^d^,x{x) = x{x)- 


(6.55) 


Going to the momentum representation x{p) before, this becomes 

iB-^Af^Pij,x(,p) = Xip)- 


(6.56) 











6.2 The Connection Between Spin and Statistics 


73 


Now consider the matrix (of finite order) iB~^= M, which must satisfy 
its algebraic minimal characteristic equation M” + + • • • + a„ = 0, 

where the a„ are numbers. From Eq. (16.201) we know that L~^ L = 

(^B~^A''). If we insert the proper combinations L~^L in each term, and 
use the relation 

L-i {B-^A^^p^) L = B-^A-if^,p^ = B-^A^p^, (6.57) 

then the result of a Lorentz transformation, giving us back the identical 
minimal equation (as it must), shows that the aj are Lorentz invariants, or 
functions of invariants. Since the only 4-vectors available are the p^, we take 
Furthermore, the minimal equation must be valid indepen¬ 
dently of the value of the numbers i.e., it must be an algebraic identity 
in ; the coefficients aj{—p^p^) must then be of the form {—p^p^Y times a 
numerical factor c^ , which is independent of the p^. We can identify the two 
possible cases, corresponding to the degree of the minimal equation being 


either even (n = 2k) or odd {n = 2k + 1); in either case the power to which 
{B~^A^^Pfj) is raised must decrease in steps of two, 

n = 2k: + (-p^p'") Ci ^-b ... 

+ {-PiiP^)^ Cfc = 0, (6.58a) 

n = 2k + l: -|-(-p^p'^) ci ^ -I-... 

+ i-Pf^P^)'" Cfc {iB~^Af^p^) = 0. (6.58b) 

Note also that the numbers Cj must be real, since 

B'f = B, (iA^^)^ = iAf^, (6.59a) 

and {iB~^A^Pfj)^ = iA^B~^Pf_i = B {iB~^A^^pY) B~^. (6.59b) 


Taking the adjoint of the minimal equation then corresponds to making a 
similarity transformation on the matrices and complex conjugating every Cj 
term. Since each B can be combined with a B~^ term, and since we must 
still have the same unique minimal equation, it follows that each cj = c*. All 
of the above is a direct inference from the requirement of Lorentz invariance. 

If we now apply (I6.58al) and (I6.58bl) to the field x(p), and use (j6.56l) we 
obtain for either case 

1 + (-Pt^P^) Cl + (-PtiP^)^ C2 H-h {-p^p^)^ Ck 


X{p) = 0. (6.60) 
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6.3 Fermi-Dirac Fields of Spin 1/2 

Consider now only the simplest case of fc = 1, which corresponds to a minimal 
equation of degree 2. Then Eq. (I6.60|) becomes 

[1 + Cl] X = 0, (6.61) 

and interpreting ci as —for a non-vanishing y we have the familiar 
relation PfiP^ + m? = 0. 

We now go to the simplest case of all, that for which n = 2 in Eq. (I6.58a|) ; 
this then becomes 

{iB-^A>^p^)^ + = 0. (6.62) 

Writing 

B = m/3, = ia^, and /3“^a^ = (6.63) 

then is Hermitian, and 7 ^^ = 13^^ Eq. (I6.62|) becomes 

h'^Pfj.f = -PtJ.p^ = -g'^'Pt^p^ = 'y^YPfj.Pi^, (6-64) 

and since only the symmetric combination enters here, we have 

= (6.65) 

as a necessary condition for Eq. (16.6211 to be satisfied algebraically by 

Before proceeding further, we remark that the construction of 5 matrices, 
satisfying (16.651) . each of a definite symmetry (to represent the 7 ^ and 7 ® = 
7 *^ 7 ^ 7 ^ 73 ), can be achieved in only one way: Three of the matrices must be 
symmetric, and the remaining two antisymmetric. This statement, which is 
easily verified, together with (I6.62p shows that we can construct matrices 
satisfying (16.651) only for Fermi-DIrac fields, as follows: Since is real, and 
/3 imaginary, i.e., all the are imaginary. From (I6.65P 

we note that 

( 7 °)^ = -500 = -^1 ( 7 '=)^ = = - 1 . ( 6 . 66 ) 

Defining = iB^ where the B^aie then real, it follows that 

(r°)^ = -i, {r^f = +i, 

E E (6.67) 

n n 

We are looking for matrices of definite symmetry; equating j and £, (I6.67|l 
can be satisfied only if B^j = — B^j = i.e., 7 ° is antisymmetric 

(and imaginary, and therefore Hermitian), and the three 7 ^. are symmetric 
(and imaginary, and therefore skew-Hermitian). 7 ® is then antisymmetric and 
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real. It then follows that {, 8 , 7 ^”} = [, 8 , 7 °] = 0, and since (3 and 7 ° have the 
same properties, we may identify them: /? = 7 °. 

If we now attempt to repeat this for Bose-Einstein fields, then the are 
anti-symmetric and imaginary, /3 is symmetric and real, all the 7 ^ are the 
same as for the Fermi-Dirac case, and again {/3, 7 ^} =0. But this is a direct 
violation of the requirements that there be but three independent symmetric 
matrices, satisfying (I 6 . 66 |) . Equation (I6.65L from which the spin ^ formalism 
is obtained, therefore refers only to Fermi-Dirac fields. 

The Lagrangian for the Fermi-Dirac spin ^ field is then 

1 777 

£ = - [' 0 , ia^d^ip] - — [ 0 , / 30 ]. ( 6 . 68 ) 

More precisely, this is the Lagrangian for the uncharged Dirac Majorana 
field—in order to represent charge we shall later have to double the num¬ 
ber of components of 0. For notation, we introduce the use of the dot • to 
symbolize the proper symmetrization brackets to be used for Fermi-Dirac or 
Bose-Einstein fields. Then 


L = 0.fa''i9^0 — TO0./30. (6.69) 

The field equations are then the familiar ones 

+ mf^ 0 = 0, (6.70) 

where, since 0 is Hermitian, the adjoint of these equations must, and do, 
yield identical equations. 

To obtain the commutation relations, we use: 

G'v' = y dafj.xA^Sx J daoiis°Stjj, (6.71) 

where s° = = i (/3 = 7 °), and form 

[X,^^] = ^^0(a:), (6.72) 

or 

J dao S-ip{x') = (6.73) 

where the factor ^ enters because the 0 are “non-canonical,” and []—>■{} 
since Sip anti-commutes with 0. This then yields 

{0(a;),0(x')} = ^ (s°) ^ S^°'>{x - xo), (6.74) 


or 
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{'ljja{x),tpf!{x')} = ^Sap6^°\x - x'), 


(6.75) 


where x and x' are points on the same space-like surface, and (s°) ^ = —i. 


6.4 Spin Zero and One 


Let us return to the characteristic equation (I6.58bl) of the matrix (iB . 

We will choose /c = 1, and examine the odd polynomial in , which 

gives the possibility of describing a particle of zero mass. This is necessary 
to describe the electromagnetic field. 

Define the four vector /3^ = imB~^A^. Then we have 

+p''Pu W^Pf,) = 0 , 

or -f g^''P'^)p^PaPv = 0. (6.76) 

From this equation, which is an identity in we can make statements 
only about the symmetric part of the matrix products. If we completely 
symmetrize the matrix factor with respect to /i, cr, and it must vanish. 
This is the sum of three terms of the form 




(6.77) 


which are the familiar Kemmer-Duffin commutation relations [Petiau 19361 
IDiiffin 19381 IKemmer 1939) which are used to describe a particle of spin zero 
and one. All the j3 matrices are singular. These matrices have 126 independent 
elements and hence are reducible to three sub-matrices of dimensionality 10, 
5, and 1; the sub-matrix of dimensionality 1 is trivially the null matrix; the 
one of dimensionality 5 and rank 2 describes a particle of spin 0; the matrix 
of dimensionality 10 and rank 6 describes a particle of spin 1. 

The Lagrangian can be written, by choosing B = m, as 


C = — m(j).(j), (6.78) 

which yields the equation of motion, 

(j) = 0. (6.79) 


Multiplying through by (3'^(3'^ jda\d^, and symmetrizing with respect to cr 
and p, using the commutation relations to reduce the triple matrix products 
and using the wave equation, we finally find 


-dj-d'' + m 
i i 


<^ = 0 , 


(6.80) 
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which verifies what has been put into the theory as the connection between 
energy and momentum. Thus each component of </> satisfies the Klein-Gordon 
equation. 

The generator is 

~ y (6.81) 

The commutation relations are developed from 






-Scj}{x), 


(6.82) 


which yields 

[^a(a;), (/3°(()(a;')) J = - x'), (6.83) 

where is a diagonal matrix having six ones and four zeros along the diago¬ 
nal, which refer to the independent and dependent components, respectively. 
This equation cannot be solved for the commutation relations because /3° 
is a singular matrix. It eliminates the dependent components of </> from the 
commutation relations. 

If we multiply the field equation by (l — (/3°)^) and use the commutation 
relations, we find 


[/3'=/3°/3°9fe + fm(l-/3°/3°)]</' = 0, (6.84) 

which is independent of time, and is the equation expressing the dependent 
components of (j) in terms of the independent components. 

Let us re-examine the spin zero and spin one fields. Integral spin fields may 
be described in terms of ordinary tensors. They do not require the special 
spinor properties of the half-integral spin fields. We have seen that the spin 
zero representation of the Kenner-DufHn algebra has a dimensionality of five. 
If the only tensor in addition to the field components which is introduced to 
form bilinear products in a scalar Lagrangian is the four-divergence, we must 
describe the field with tensors differing in rank by one. We shall construct the 
spin zero field with a scalar and four-vector as the necessary 5 components, 
and the spin one field with a four-vector and an anti-symmetrical second rank 
tensor as representing the 10 components. 


6.5 Spin Zero 

The Lagrangian for a spin zero field is 

1 777 

^=2 - y <(''') ■ (6-85) 


The field equations are determined by varying (j) and (j)^, 
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= mcj), d^<j) = ( 6 . 86 ) 

which imply 

d^dncj) = md^(j)n = d^dncj)’' = —d'^d^di^cj) = (6.87) 

m 

again yielding the Klein-Gordon equation. 

The generator is 

G = J daf^^{(j)S(l)'^ - (j)>^S(j)) ^ J dao^ ■ (6.88) 

Thus only (jj and (j)^ are independent variables. This is reflected in the fact 
that the Klein-Gordon equation is second order, and hence we must specify 
both the wave function and its time derivative. 

The equations of motion have to be examined to see if the other com¬ 
ponents of the field are determined. We have from the field equations, as 
equations of motion. 


do4' = m(i)Q d^cjP = mcj} — dkfj)^, (6.89) 

and the following which is not an equation of motion, but does show that the 
(f>^ are determined in terms of the two independent components, 

dk4> = 'mcj)k. (6.90) 

From the generator we see that we have one set of canonically conjugate 
variables, cj) and (f>Q. This means that the field has only one internal degree 
of freedom, and must be a spin zero field. This pair of conjugate variables 
must obey the canonical commutation relations (at equal times), as is easily 
verified by using the generators for their respective change 

[(j){x),(l){x')]= [(j)°{x),(j)°{x')] =0, [(j)°{x),(l){x')] = i6^°\x - x'), (6.91) 

where — x') means (5(r,r') on uo- The commutation relations obeyed 

by the dependent field components are derived from the field equations and 
these; for instance 


- [(I)°{x),(j)k{x')] = -—dhS^^^x - x'). (6.92) 

i m 

These commutation relations can be written by inspection in a form not 
referring to any particular coordinate system (where now S^^\x — x') means 
(5(r, r') on a^), on the spacelike surface cr^. 
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[(l){x),4>{x'))] 


0 ^[(j)^{x),(l)ix')] = S^^^x — x'), 

_1 \df^S^'''>{x-x')+d''S<-^\x-x') 


(6.93) 


6.6 Spin One 

The description of a spin one field requires the use of a ten-component wave 
function. As we shall see, not all of these components are independent. We 
shall use a four-vector and an anti-symmetrical tensor of rank two. The La- 
grangian is chosen to be 

£ = i - I . (6.94) 

The factor of 1/2 in the last term is inserted because the unrestricted sum 
over fj, and v counts each component twice. The field equations are 

dfj,4>v - d^(f)^ = mG^y, dfj.G^'-' = mcj)'^. (6.95) 

The generator is 

G = J dao^ (</>fc5G°'= - G^’^Scj.k) . (6.96) 

Thus only six of the ten components can be varied independently, <j)k and 
G°^. Let us see if they have equations of motion and if the other 4 components 
are determined in terms of the independent ones. The equations of motion 
are 

do4>k = + mGofe, 5oG°'= = m(t>^ - diG^\ (6.97) 

for the independent field components. Also 

dk4>t - dt(j)k = rnGki, mcjp = dkG^° = -dkG°'‘, (6.98) 

which determine the four dependent field components. Thus all ten compo¬ 
nents of the field are determined. 

From the generator we see that we have 3 sets of canonically conjugate field 
variables; thus the field has 3 internal degrees of freedom. This corresponds 
to the three sub-states of a spin one field. 

The commutation relations obeyed by the field components are obviously 
the canonically conjugate relations, 

[</.'= (x),/(x')] = [G°'=(x),G°^(x')] =0, 

-\G°'^{x),Mx')] = S^S‘'°\x-x'). (6.99) 

i 
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The commutation relations obeyed by the dependent components can be 
realized from the above, by using their definitions in terms of the independent 
components. For instance, 

- \(jp{x),cj)^{x')\ =- S^dkS'^^^x — x') = - diS^°\x — x'). (6.100) 

i m m 

These relations can be generalized to refer to an arbitrary coordinate sys¬ 
tem. 


-[ck^{x),r{x')] 

I 

I 

i \Gf^'^{x),G^^ix')] 
i ■' 


- di^S^''\x-x') + d‘'S^>^\x-x') , 

ml J 

{x - x') - (a; - x'), 

(d''5‘'^\x - x') + d^5‘'''\x - x')) } 
plus antisymmetrical terms .... (6.101) 


6.7 Electromagnetic Field 

The electromagnetic field is a spin one, massless field. A re-definition of the 
field variables in the limit as the mass approaches zero will be made in the 
spin-one Lagrangian. Set 

( 6 . 102 ) 

VTO 

Then the Lagrangian ()6.94p becomes 

>^ = \ + \f^''F^,. (6.103) 

The field equations are 

= 0, (6.104) 

and the generator becomes, in the local coordinate system, 

G = J da^{AkSF°’^-F°’^6Ak). (6.105) 

Thus, as in the case of a spin one non-zero mass field, only six of the field 
components can be varied independently. Their equations of motion are 

doF^’^ = -diF^^. 


doAk = dkA[)+ Fofc, 


(6.106) 
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We must now examine the rest of the field equations to see if the dependent 
components are determined, 

dkAi - diAk = Fu, dkF°^ = 0. (6.107) 

Thus F^’- is determined, but because of the vanishing of the photon mass, 
Aq is not determined and a further restriction is placed on the electric field, 
F^k _ This is exactly the freedom of making a gauge transformation, which is 
allowed by the vanishing of the photon mass. That is, 

Ajj^ —> Ap^ + dpA. (6.108) 

It is convenient to split the vector part of Ap into a longitudinal part, which 
is the gradient of a scalar, and a transverse part, which is the curl of a vector. 
Thus we see that the gauge transformation affects only the longitudinal part 
of Ak- F'^^ is purely transverse because it has no divergence. 

The equations of motion also break up into longitudinal and transverse 
components, 

doA\ = dkAo doAj. = F^k- (6.109) 

The first equation is entirely consistent with the possibility of making a gauge 
transformation which implies 

A\. = dkA, Ao = doA. (6.110) 

By setting 

Ak = dkA + Aj ( 6 . 111 ) 

in the generators and integrating by parts, or equivalently adding a surface 
term to the Lagrangian, we see that only the transverse components of Ak 
appear in the generator, 

G = J dao^ {AJSF°'^ - F°HAl) . (6.112) 

Since F^^ and A'^ appear as canonical variables, we can impose the canon¬ 
ical commutation relations 

[Al{x),Ajix')] = [F°>^ix),F<^^{x')] = 0 , (6.113a) 

[F^'^{x),AJ{x')] = — x')^ , (6.113b) 

where only the transverse component of the term on the right side of (I6.113bl) 
can appear. If we define a term that subtracts off the longitudinal part of the 
delta function, the commutation relation can be written as 

[F°'=(x), Aj(x')] = S^S^°\x -x')- d’^d'eV°{x - x'). 


(6.114) 
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The divergence of this purely transverse equation must vanish, hence 

SgdkS^'^\x — x') = —d^dkd(jf‘{x — x'), (6.115) 

which can be satisfied if 

V^V'^{x-x') = -5^'^\x-x), (6.116) 


or 


V^{x-x') = 


1 


(6.117) 


47 r|x — x'l ’ 

which shows that the longitudinal part of is intimately connected with 
the Coulomb potential. 


6.8 Introduction of Charge 

It was previously mentioned that the description of charge requires a doubling 
of the number of component fields x- 

The simplest case to consider is the case of one internal degree of freedom. 
Calling the two basic fields X(i) and X{ 2 )^ -C [Eq. (16.81) ] becomes 

^ = X(iyA^d^,X(i) + X{2)-A^d^X(2) - 'H(X( 1 ), X( 2 )). (6.118) 

Since both X{i) and X( 2 ) have identical space-time properties, the kinemat- 
ical portion of £ is invariant under rotations and reflections in a new two- 
dimensional space composed of the components X(i) and X( 2 ) in accordance 
with the sum-of-squares notation of (16.1181) . The rotations in this space, given 
by 


X( 1 ) = cosAx(i) -ksinAx(2), 


X( 2 ) =-sinAx(i)+ cosAx( 2 ), (6.119) 


may be concisely expressed as 


X = 




X, 


( 6 . 120 ) 


where 



( 6 . 121 ) 


The basic improper transformation (reflections) can be taken either as 


X(i) — X( 2 ), 
or 

X(l) = X( 2 ), 


II 

IX 

( 6 . 122 ) 

X( 2 ) — “X(l)- 

(6.123) 
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Re-defining as 

(T 1 ) ■ 

allows us to write (|6.118|) in matrix notation 

C = x-A^d^X-nx). (6.125) 

where we have assumed that % has the required invariance under reflections 
and rotations in this space. Use of (16.1201) and = —q is then sufficient to 
guarantee the invariance of C under rotations. The reflection transformations 
in (16.1221) can be compactly written as 

X = C'e*5«x X = Cx, (6.126) 

where 

C = (; J) (6.127) 

Then the relations = C~^ = C show that C is invariant under reflections. 

The proper transformation (16.1201) can be evolved from a sequence of in¬ 
finitesimal transformations of the form 

X = (X+i^M)x = X-^X Sx = -'iSXqx, (6.128) 

and incorporating this infinitesimal variation into the action principle will 
yield a quantity which is conserved. The notation will suggest the interpre¬ 
tation of this quantity as electrical charge, but until the manner in which it 
appears in interactions between the various fields is specified, the nature of 
this charge is irrelevant. 

If we now imagine that <5A of |6.128| is a continuous function of space time, 
the finite change in A which builds up between an initial t{a 2 ) and final 
t{ai) times corresponds to different successive choices of x The stipulated 
invariance of C then implies that 5\C = 0. If (5A were constant, 5\C would 
vanish trivially, in accordance with (16.1201) and (I6.126|) ; we obtain something 
new since the dependence of 5\ on introduces terms depending on 9^^A: 

5\C = -ix-Af^qxdfj,5X = j^9^(5A = [j^(5A] - SXdfj,j^, (6.129) 

where 

= -iX-A^qx = -ixA^qx = Xi2)A^Xii) “ X(i)^''X(2), (6.130) 

The action principle now requires that the coefficient of 6X vanish, yielding 
the differential conservation law 


= 0 . 


( 6 . 131 ) 
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The term 9^ gives the generator for changes on the boundary space-like 

surfaces CTi. a 2 - Taking 5\ as constant over each surface, which corresponds 
to different, but definite, choices of the x o’i,2i tve have 

Gx=5\ J da^f = QSX. (6.132) 

Further, if (5 A((Ti) = ( 5A((T2) then it is obvious that we get the same description 
as we would have obtained with (5A = 0 in both surfaces, i.e., 

[Q(o' 2) — Q(o'i)] ^A = 0, and therefore Q(cri) = Q{(y 2 )- (6.133) 

This result also follows from (j6.13111 : the quantity Q, called the “charge”, is 
therefore conserved. 

Obtaining the commutation relations, from (I6.128|) and (16.1321) we have 

^ [X, ^] = 7 [X> Q^A] = -iqxSX, (6.134) 

or 

[x,Q] = 9X, (6.135a) 

or 

[X(1),Q] =-*X(2 ), [x( 2),Q] =-*X(i)- (6.135b) 

The generalization of (|6.134l) for a finite rotation is 

which is analogous to translation induced by the momentum operator 

[x,Piy] = 

e-*“'^x(x)e“-^ = = x(^ + «)• (6-137) 

There must also exist a unitary operator, in analogy with ()6.136|) . which 
accomplishes the improper transformations of equation (16.1231) : 

C-^xC = cx, or . (6.138) 

(C "X(2)C = -X(2) 

From (I6.132P and (16.13011 we then have 

C-igC = -Q. (6.139) 


Thus, C has the interpretation of a charge reflection operator. Since = -1-1, 
its eigenvalues are C = ±1. From (|6.13911 we see that [Q^C] yf 0 and the 
two operators cannot be simultaneously diagonalized, (in the representation 
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of X(i):X(2):C is diagonal) except for states of zero charge. If we arbitrarily 
assign the vacuum state the eigenvalue C' = +1 , then this state of zero 
charge, Q' = 0, has both a definite charge and charge symmetry. 

To obtain a state (other than Q' = 0 ) of definite charge, where the 
operator Q is diagonalized, we must utilize non-Hermitian operators. Define 

X+= X(i) - *X(2), X-= X(i)+*X(2), (6.140) 

and then a simple calculation shows that 

[x+,Q]=X+, [X-,Q] = -X-, (6.141) 

indicating that in this representation Q is diagonal. However, C is no longer 
diagonal: 

C-\±C = Xt- ( 6 - 142 ) 

From (16.1411) we see that if we have a state IQ') representing a definite 
charge, then 


q(x+|Q')) = (Q'-i)(x+IQ')) 

Q(x-IQ')) = (Q' + I)(x-IQ')) (6.143) 


showing that x+ creates a state for which the eigenvalue of Q is Q' — 1, and 
conversely for x- These, then, represent charge annihilation and creation 
operators respectively. If we imagine that every physical state can be created 
by the action of these operators on the vacuum |Q' = 0) then the only possible 
values of Q' are 0, ±1, ±2, ... etc. (Note that charge annihilation represents 
the destruction of positive charge or the creation of negative charge, and vice 
versa.) To connect this with our previous particle description in Sec. 16.11 
where we had x^"*"^ X^~\ obviously have the four possibilities: 

• X^"^ destroys a particle and decreases charge by 1, 

• X- ^ creates a particle and increases charge by 1. 


These are inverse operations: ^X- = X^"^- 

• X^^ destroys a particle but increases charge by 1. 

• X^^ creates a particle but decreases charge by 1. 


These are also inverse operations: (x^^) = X+ This exhaustive descrip¬ 
tion permits us to describe particles that carry charge. 

In terms of the non-Hermitian variables defined by (16.1401) . the Lagrangian 
(|6.118|) becomes 


^ = \ [x--A^d^x+ + X+-A^df,x-] - 


( 6 . 144 ) 
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and the generator is 

G = y [x-A^^Sx+ + x+A^Sx-], (6.145) 

which shows that x+ and X- are a canonically conjugate set of field variables. 
The current (16.1301) is 

\ [X+^^'X- - X-^'^X+l ■ (6.146) 

For variables of the second kind the commutation rules are 
{x+ix),X-{x')} A° = i6^°\x - x'), 

{X±(a;), X±(a;')} = 0- (6.147) 


6.9 Quantum Electrodynamics 

As an illustration of the use of non-Hermitian fields, we now consider the ordi¬ 
nary charged Dirac field. We start with and '0(2) as the two 4-component 
Hermitian fields which form the simplest charged field. Let 

0(1) - #(2) = 0+ = 0, (6.148a) 

0(1) -|-10(2) = 0- = 0^- (6.148b) 


There is an artificial asymmetry here depending on what is defined as 0 and 
01. The commutation rules are now 


|0a(a;),0^(x')| = Sai3S^°\x - x'), 
{0a(x),0;3(ai')} = 0 = |0),(a:),0^(ai')} • 

(6.149) 

Since in this case, = i, the Lagrangian becomes: 


£ = - ((?p0) — i (i9ii0^) q;'^ 0] — m0)’/30. 

(6.150) 

The current is then 

= 0fa^0. 

(6.151) 

Setting 

, ilA P = 0, 

(6.152) 

the Lagrangian is 


^ ^ [0.7^fcli,0 - icli,0.7'^0] - m0.0. 

(6.153) 




6.9 Quantum Electrodynamics 


87 


If we wish to interpret as the electric current, we can write down the 
Lagrangian for the electromagnetic field plus the charged Dirac field and 
specify the coupling so that the usual field equations arise, in which acts 
as the source for the vector potential. The form of the coupling term is of 
course limited by Lorentz invariance and the spin i and spin 1 algebra. Thus 
the most general Lagrangian which can be formed in this case is 





(6.154) 


The symmetrization between and is necessary because in general, the 
sources for the are partly the and hence these need not commute. The 
last term, a Pauli moment, even though it appears covariant, may, in fact, 
not be covariant due to its operator properties. These show up only in higher 
order terms in the perturbation expansion. 

The Lagrangian given by equation (I6.154p is invariant under the rotation 
and reflection of X(i) X{ 2 )^ which correspond to the replacements 


A^^ A^ + ^d^A{x), (6.155) 


and to the interchange of '(/' and 'ip, respectively. Note that in the second 
case—that of reflection—we must also replace —Afj,, F^i, —>• —F^y, in 

order to have C unaltered. 

We obtain the equations of motion and generators by variations 



(6.156a) 


Sip : Ip [7^ (i9^ 

SF^y : d^Ay - dy/ 

6A^ : a.F'^" = ei 


V' [ 7 ^ (*9j - eAf,.) 


eA^.) + m] = 0, 


(6.156b) 

(6.156c) 

(6.156d) 


^^Ay dyApi Fpj^y , 

dyF>^'' = ePp.x^iP = f 


The generators are 



(6.157) 


and by adding appropriate 4-divergences, this becomes 



(6.158) 
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in a local coordinate system. In this form, the independent variables appear 
to be only ip and the . ip is definitely an independent variable; its equation 
of motion is Eq. (|6.156al) . To re-examine the electromagnetic field, we write 
Eqs. (I6.156cl) and (I6.156dl) as 

doAk = dkAo + Fok, Fm = dkAi — diAk, 

g^F>co = f - deF^°, dkF°’^ = f. (6.159) 

Since we have the freedom of making a gauge transformation, wherein the lon¬ 
gitudinal components of A>^ are then arbitrary, it is advantageous to re-write 
equations (16.15911 in terms of transverse (T) and longitudinal (L) components: 

(1) = (2) doF’^^^ = - d,F^\ 

(3) Fke = F^pf^ = dkAi — dgAk, (4) d^A^ = dkA^ F Eq),, 

(5) doF’^°^ (6) dkF°'^^=f. (6.160) 

Here, F^f. ^ 0 since the electric field is no longer divergence-free. The solution 
of item (6) is immediately 

p 0 kL^_Qk J da^°^V(^o)ix - x')j°ix'), (6.161) 

demonstrating that the longitudinal part of the electromagnetic field is not a 
kinematically independent quantity, but depends on the of all the charged 
particles present. Note that item (5) is not an independent statement, since 
it is the result of the conservation law for j^{x). Since A\ = dkA{x), using 
item (5) shows that 

Aq = doA{x) + J {x — x')jo{x'), (6.162) 

i.e., Aq is completely arbitrary, corresponding to the freedom of choice of 

Mx)- 

The generator in equation (16.1581) can now be written as 

G = J da [pi-f°Sip - F^'^'^SAl - F°’^^SAl] . (6.163) 

By the addition of a surface term, the last term in the integrand in (I6.163|) 
becomes equal to = j^SA, which does not refer to the kinematically indepen¬ 
dent part of the electromagnetic field, and gives the generator for the change 
in the Dirac field when an infinitesimal gauge transformation is made—See 
Eq. (j6.132|l . The electromagnetic field commutation rules are again as given 
in equations (I6.113al) and (I6.113bl) : the anti-commutation relation for the 
Dirac fields works out to be 

{ip{x),ip{x')} = 'y^S^^\x—x') = —jo6^^\x—x') —>■ —"f^5^^\x—x). (6.164) 
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Using the non-Hermitian ip, ip, the effect of the charge reflection operator 
C is to interchange the fields: 

C-Vc = ^^ C-V^C=t/>. (6.165) 

In order to maintain the invariance of C, we must now also require that this 
unitary operator reverse the sign of the 

C-^Af^C = -A'". (6.166) 

Again, = +1, and therefore its eigenvalues are C = ±1. From the argu¬ 
ments of the proceeding section, use of the transformation: 

C-iQC = -Q (6.167) 

permits a classification, in terms of the quantum number c', for systems of 
zero net charge. For example, consider the production of photons by repeated 
application of the operator A (actually A^~\ a photon creation operator) to 
the vacuum state |0), characterized by c' = -1-1, and assume that no coupling 
terms are present in C (i.e., e = 0). If A is applied n times to |0) then by 

C (A • • • A) |0) = (-1)” (A • • • A) |0), (6.168) 

and we have a state |(A' • • -A')) whose eigenvalue of C depends on the (even 
or odd) number of photons present: C = (—1)" = ±1. Whichever -|-1 or 
-1 eigenvalue of c' we begin with persists as the coupling is turned on, even 
though what is now called the state of n photons is a superposition of many 
different states. 

As an example of this, let us consider the decay of positronium. In the 
absence of the coupling term between A^ and ip (i.e., e = 0), the system of 
(e+ + e~) is stable, and has a definite eigenvalue c'. When the coupling is in¬ 
troduced, the only permitted final states will be those of the same eigenvalue. 
From our previous interpretation, ip is the operator which creates a particle 
e~ and annihilates , and conversely for ip’’^. To obtain a state with e~ at x 
and e’*' at x', we form (ip^{x).ip{x')) |0), and recognize that the actual state 
is a superposition of these with the amplitude wave function ip{x,x'). Using 
(16.1651) we then have 


C|e+e“) = C J J {dx){dx')ip{x,x') {ip\x).ip{x')) |0) 

= —J J {dx){dx')ip{x,x') [ip\x').ip{x)) \0) 

= —J J {dx){dx')ip{xPx) [ip^{x).ip{x')) \0). (6.169) 


If the wave function ip{x, x') is symmetric, then d = —1; if ip{x, x') is anti¬ 
symmetric, d = -|-1. As the coupling is now turned on, this no longer remains 
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the exact state, but the d value remains the same. Thus we have the selection 
rule for positronium: from a state it can decay only into an even number 
of photons (c' = +1, tjj{x,x') is antisymmetric), from a state it can decay 
only in to an odd number of photons (c' = —1, il){x,x') is symmetric). 


Chapter 7 

Nonrelativistic Source Theory 


The following is based on lectures given at Schwinger by Harvard in Spring 
1969, as transcribed by the author. The goal was to construct a general theory 
of particles, in a nonrelativistic context. As such, this provides a transition 
between nonrelativistic quantum mechanics and source theory, the general 
development of which was given later in Schwinger’s three-volume treatise 
[Schwinger 1970b| [Schwinger 1973[ [Schwinger 1989a| . 

The measurement symbol, or projection operator, which forms the basis 
for Schwinger’s approach to quantum mechanics [Schwinger 2001[, 



(7.1) 


represents a idealized process in which the state b' is annihilated and a new 
state a' is produced. Wouldn’t it be useful to separate these processes? Recall 
how oscillator states, as discussed in Chap. [H were created by forces. We 
generalize to sources which can create and destroy particles. A scattering 
process, which occurs in a more or less localized scattering region, can be 
abstracted into a two-stage process, in which first an incoming particle is 
absorbed, and then a new particle is created, as sketched in Fig. 17.11 In the 
individual processes the source K acts to create or annihilate the particle. 

What does this mean quantum mechanically? We must describe the pro¬ 
cesses by probability amplitudes; we want to construct everything from the 
one particle production mechanism, 


(lp+|0_)^, (0+|lp-)^, 


(7.2) 


where the first amplitude represents the process in which a single particle 
state of momentum p is created, where before the source acts only the vacuum 
state is present, while the second represents the process in which a single 
particle of momentum p is absorbed after which the vacuum state is present. 

Because the processes occur in space and time, the source must be a func¬ 
tion K (r, t ), which exhibits a certain degree of localizability. How does the 
effectiveness of the source vary with different degrees of freedom? The comple- 
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Fig. 7.1 A scattering process in which a particle of momentum p 2 scatters into a 
state of a particle with momentum pi can be thought of as a combination of two 
processes: One in which the particle is created by a source K in the state pi and a 
second in which the particle of momentum p 2 is absorbed by the source K. 


mentary measure is the corresponding function in momentum space K{p, E), 
where, nonrelativistically, E = /{2m). We expect the relationship between 

the production amplitude and the source function to be, at least for a weak 
source, 

(lp+|0_)^ = y^^H)i^(p). (7.3) 

Here (dp) = dpidp 2 dp 3 . The —i factor is purely conventional for later conve¬ 
nience. The square root of the momentum-space element is present to prop¬ 
erly account for the density of states in the continuum picture. To compute 
the annihilation amplitude (0+|lp_), we can use orthogonality, (lp_|0_) = 0. 
This must be maintained by the dynamics. The completeness relation 

1 = |0+)(0+l + X/ l^p+)(lp+l + X/ |lp+lp'-f)(lp+7p'-fI + • ■ •; (7.4) 

p p.p' 

implies 

0 = (lp_|0+)^(0+|0^)^ + ^(lp_| V+)(V+|0_)^ + ..., (7.5) 

p' 

where we shall consider a weak source, so we will drop the higher terms. To 
lowest order (in powers of the source) 


(0+|0_)° = l, (lp-|V+)° = dpp,, (7.6) 


SO 


o = (v|o+)^ + (ip+|o_)^ 


( 7 . 7 ) 
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or 

(0+|lp_)^ = -[(lp+|0_)^]*. (7.8) 

Therefore, the effectiveness of producing or absorbing a particle by a weak 
source is 


{lp+\0-)^ = (0+|lp_)^ = ySl(_*)K(p)*. (7.9) 

These equations, in fact, define what we mean by a source. 

Now we need to seek the relation to the space-time description. What 
happens when a source is displaced, 

K{r,t) = K{r + R,t + T). (7.10) 

Relativity (here Galilean) means that the same effect occurs by displacing 
the space-time coordinate system to which the initial and final states are 
referred, the generator of such an infinitesimal displacement being 


G = P ■ 6e- HSt, (7.11) 

which implies that the displacement operator is 

U = (7.12) 

Thus the 1 particle states change according to 

(lp|^(lp|e =e*^ (Ipl, (7.13) 

while the vacuum state is unchanged, 

(0+1 ^ (0+|. (7.14) 

Thus if the one particle amplitude is proportional to K{p), 

(lp+|0_)^~K(p), (7.15) 

that of the displaced source is 

(lp+|0_)^ ^ e'P'^-*^^if(p) ^ K{p). (7.16) 

This implies that K{p) is obtained by Fourier transforming K{r,t), 

K{p) = K{p,E) = j {dr)dte-^P"^+^^^K{r,t). (7.17) 


In general, one wants to remove the connection between E and p. 
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What about strong sources? Remember for the oscillator, the most basic 
object was the ground-state persistence amplitude (0+|0_)^. The latter con¬ 
tains the process of the exchange of particles between temporally separated 
sources, as illustrated in Fig. 17.21 Here a single particle is emitted by source 

- 0 + 






K2 


0 - 


Fig. 7.2 Exchange of a particle between a source K 2 and a source Ki, where the 
latter is localized later than the former. Before the first source acts, the system is in 
the vacuum state, as it is after the second source acts. Time is imagined as plotted 
vertically in these “causal” diagrams. 


K 2 and absorbed later by the source Ki. Before and after either source acts, 
the system is in the vacuum, no-particle, state. We imagine the total source 
to be the sum of the two components, 

K{x) = Ki{x) + K 2 {x). (7.18) 


The decomposition shown is called a causal arrangement. Because we are 
so far considering weak sources, when we insert a complete set of states at 
an intermediate time between the action of the two sources, the vacuum 
persistence amplitude is 


(0+|0_)^ = (0+|0_)^^(0+|0_)^^ +^(0+|lp_)^^(lp+|0_)^= +.... (7.19) 

P 


Now using the one-particle creation and annihilation amplitudes (Tza , to¬ 
gether with the Fourier transform (17.171) we see that the one-particle exchange 
term here is 


— i 


{dr)dt{dr')dt'Ki (r, t) 


—i 


(^P) iv(r-r')-iE{t-t'') 

(27r)3 


K2{Y',t'). (7.20) 
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But the source is a unitary whole—results can depend only on the total 
source K and not its parts. This is a statement of the uniformity of space 
and time. This will introduce terms that will refer to each component source 
separately; what we don’t want is a term that involves and Ki] therefore, 
we define the retarded Green’s function 

= j ^ (7 2I) 

in terms of which we infer 

(0+|0_)^ = l-i j{dr)dt{dr')dt'K*{r,t)G{r-r',t-t')K{r',t'). (7.22) 

The Green’s function satisfies the differential equation 

G(r — r\t — t') = S{t — t')6{r — r'). (7.23) 

Now we want to remove the restriction to weak sources. Suppose we have a 
beam of noninteracting particles, detected for example by spatially separated 
sources, as illustrated in Fig. 17.31 Each particle is produced and detected by 


0 + 



- 0_ 

Fig. 7.3 Exchange of noninteracting particles between spatially nonoverlapping 
sources. Particles emitted by the lower source of ifGl are only detected by the upper 
source of the same set. It is arranged that there is no cross-coupling. 


(jv) 

dt 2m 


a single pair of weak sources K^°‘\ a = 1, 2, 3,.... There is no interaction 
between different pairs of sources. Thus the vacuum persistence amplitude 
for this arrangement is 
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\ f 1 

(O+lO-)*” = 1 — i j {dv)dt{dv')dt'K*{v,t)G{v— y' — t')K{v',t') 

(7.24) 

because each source is weak. But only K, not should enter; physics 

shouldn’t depend on the channel. So 


where the last step depends upon the arrangement that prohibits cross cou¬ 
pling between the component sources. This looks just like the structure we 
saw for the harmonic oscillator (j4.91[) , except now the integrals are over space 
as well as time: 


(0+|0_)^ = ^-^IGT)dt{dr')dt'K*{r,t)G{T-r',t-t')K{T',t')^ (7.26) 


We will apply this to generalized (not weak) emission and absorption pro¬ 
cesses, but still assuming that the particles are not interacting. 


0 + 



Ki 



K2 


0_ 


Fig. 7.4 Exchange of noninteracting particles between causally separated sources, 
K 2 and Kx- 


Figure [7(3] shows the exchange of noninteracting particles between causally 
separated sources, K = Ki-\- K 2 . The total vacuum persistence amplitude is 

(0+|0_)*^ = (0+|0_)^i(0+|0_)^’^e"*-^^i'^^^ (7.27) 


The fact that there is a causal relation between the two sources means that 


G(r — — t') 


jdp) ir}.(r-T')-iE(t-t') 

(27r)3 


(7.28) 
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Define the discrete specification of the source in momentum space as 

K, = ( 7 - 29 ) 

Then the final exponential term in Eq. (17.271) is 

(7.30) 

On the other hand, if we insert a complete set of multiparticle states at an 
intermediate time, we have 

(0+|0_)^ = ^(0+|M)^^(M|0_)^^ (7.31) 

{n} 

Compare this with the expansion of Eq. (17.271) : 


( 0 +| 0 -)^ 


°° (—’iK* ^"p (—iK* 

(o+io-)^^TT y: ^ X ^ M (0+10-)^-, 


p np=0 




(7.32) 


where the occupation numbers in a given momentum cell are given by 
{n} = {u-p}. From this we infer the probability amplitudes for producing 
and absorbing particles by a strong source: 


(0+IW)^ 


n 

p 

n 



(0+|0_)^, 



(0h-|0_)^. 


(7.33a) 

(7.33b) 


As a check of this, we verify that the total probability must be unity: 

1 = ^p({n},0)^ = |(0+|0_)^pe^p (7.34) 

{ri} 


Independently, 

|(0+|0_)^P=exp 


- J{dr)dt{r')dt'K*{r, t)iG{{r — r', t — t')K{r', t') 
{dr)dt{r')dt'K*{r, t)[iG{{v' — r, t' — t)]* K{y' , t') 


(7.35) 


But the combinations of Green’s functions appearing here is 
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iG{{r-v',t-t') + [iG{{r'-r,t'-t)Y = J (7.36) 

which is a solution of the homogeneous equation, so indeed 

|(0+|0_)^P = (7.37) 

so the probability condition (17.341) is satisfied. 

The description so far of the exchange of particles between sources is a 
sort of action at a distance picture. We are often concerned with excitations 
produced by a source—a more local description. A test source is used to 
measure effects. So let us add an additional infinitesimal source, 

K{r,t) ^ K{r,t) +SK{r,t), (7.38) 

which results in the following infinitesimal change in the action, 

5W = — J {dr)dt[SK*{r,t)ilj{r,t) + SK{r,t)ip*{r,t)], (7.39) 

which defines new objects which refer to the pre-existing situation. Here 

'(/'(r, t) = J {dr')dt'G{r — r',t — t')K{r',t) (7.40) 

'0*(r, t) = J {dr')dt'K*{r',t)G{r' — r,t' — t)K{r',t'); (7.41) 

the latter is not the complex conjugate of '0, because 

G(r t'Y Y G(r' -r,t' - t), (7.42) 

G being the retarded Green’s function. The differential equation satisfied by 
the Green’s function (I7.23|l is 

G{r - = S(r - r')d(t - t'), (7.43) 

where T = —V^/(2m) is the kinetic energy differential operator. Therefore, 
the field ijj satisfies 

(^i-^-T^tPY,t)=K{r,t). (7.44) 

This is analogous to the equation satisfied by the harmonic oscillator variable 
(I4.65p . or 

- w) y(t) = Ar(t). (7.45) 


Now because 
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G(r' -r,t' -t) =S(r- r')S(t - t'), (7.46) 

the field tp* satisfies 

(-i-^^-T^r{r,t)=K*{r,t), (7.47) 

which is the complex conjugate equation. The boundary conditions are dif¬ 
ferent in the two cases: is a retarded solution, while tp* is an advanced 

solution. 

What does W have to do with action? Let us write the alternative forms 


IT = — y {dr)dt K* (r, t)ip{r, t) = — J{dr)dt tp* (r, t)K (r, t) 

= - I{dr)dtip*{r,t) (7-48) 


Combining these forms appropriately, we can write 


d 


W = j {dr)dt 

Think of this last as a functional of K, K*, ip, and ip*, so 


{i-^-T]iP- K*iP - iP*K 


(7.49) 


6W 


{dv)dt [—5K*ip — ip*6K] + 


(7.50) 


But the definition (17.401) of ip and ip* shows that the last variation is zero, 
which is a statement of the stationary action principle: 


= J {dv)dt6ip* 

+ [ {dr)dtSip 


^^^-T]iP-K 


-i--T]r-K* 


= 0 . 


(7.51) 


That is, the stationary action principle, that W is unchanged under infinites¬ 
imal field variations, supplies the equations of motion (17.441) and (17.471) . 

Fundamental to physics is the notion of the uniformity of space and time, 
that the laws of physics are independent of the locale. This is reflected in 
the indistinguishability of identical particles, which, in turn, is reflected in 
the probability amplitude of a source producing a multi-particle distribution, 
Eq. (j7.33a|) . which says all that is possible. Experimentally, we know of two 
kinds of statistics. Here rip is unlimited, so this cannot refer to Fermi-Dirac 
statistics. In the Bose-Einstein case we have stimulated emission. Let us see 
this. 
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Figure 17.51 shows the interchange of noninteracting particles between 
sources K 2 and A"i, but now with a weak source Kq in between. The to- 


0 + 



Ki 


Ko 



K2 


0_ 


Fig. 7.5 Effect of a weak source Kq acting at an intermediate time between strong 
sources K 2 and Ki. 


tal source is composed of three causally separated pieces, 

K = Ki+Ko + K2. (7.52) 

Using this causal arrangement, the vacuum persistence amplitude is 

= (7 53 ) 

Here, because the disturbance by Kq is regarded as weak, we approximate 
(0+|0_)^“ « 1, and the exponential term is expanded to first order in Kq: 

{n} 

X |l + ^[(-iiFrp)(-*i^op) + (-*iFo*p)(-*^2p)]| .(7.54) 
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Compare this with the multi-particle exchange description between the three 
sources, 

(0+|0_)^= ^ (0+|M_)^^(W+|W'_)^»({n}V|0_)^^ (7.55) 

{n}{n}' 

Now recall the connection between Kp and ({n}|0_)^, Eq. (I7.33b|) . so we 
infer 

(0+|{n})^n-*^i*p) = v^^^;TT(0+|{n} + lp)^b (7.56) 

Therefore, for a weak source, 

{{n + lp}\{n})^ = (7.57) 

The corresponding probability of creating 1 more particle with momentum p 
is 

|({n+lp}|M)^P = |XpP(np + l). (7.58) 

In the last factor, the 1 represents spontaneous emission, and the rip is the 
enhancement effect of stimulated emission. 

The corresponding analysis in the absorption case gives 

{{n}\{n + lp})^ = (7.59) 

or 

{{n - lp}|{n})^ = (7.60) 

that is, the probability of absorbing one particle is proportional to the incident 
intensity. 

Let us come back to space and time. Adopting a more telegraphic notation, 
we can write the vacuum persistence amplitude as 

(0+|0_)^ = ^ ^ {1)0(1 - l')K{l') 

-iy dldl'd2d2'A:*(l)A'*(2)G(l-l')G(2-2')A:(T)A:(2') + .... 

(7.61) 

Here the numbers represent space-time points, 1 = ri,fi, etc. Because the 
product of sources is symmetrical, we can replace 

G(1 - l')G(2 - 2') ^ i[G(l - l')G(2 - 2') -b G(1 - 2')G(2 - T)], (7.62) 

Diagrammatically, the two terms can be represented as in Fig. 17.61 This builds 
in the symmetry in the labels. The third term in Eq. (17.611) is 

-^Jdld2K*{l)K*{2)^{l)^{2), V'(1)V'(2)=V'(1,2) = V'(2,1), (7.63) 
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Fig. 7.6 Exchange of two particles between spatially and temporally separated 
sources. Bose-Einstein symmetry implies that the particles are exchanged between 
either spatially separated source. No interaction is to be inferred where the lines 
cross. 


which exhibits another characteristic of Bose-Einstein statistics. The field 
Ip, the generalization of the wavefunction, must be totally symmetric under 
interchange of the particles. 


7.1 Interactions 

Real particles interact with each other. Thus we should have processes such 
as sketched in Fig. 17.71 The particle emitted by the sources are scattered 
by interactions represented by the black box labeled T. Let us begin with a 
simpler situation, scattering from a fixed center, for example, a heavy nucleus. 
This process can be represented with only two sources, as shown in Fig. 17.81 
This diagram schematically represents the tracks of particles as seen in a 
detector. In the absence of interactions 

IF = - y K*G°K, (7.64) 

where now we have used the superscript 0 to designate the free particle 
propagator or Green’s function previously denoted simply by G. This process 
in represented by the dotted line in Fig. 17.81 Now we want to add something 
to this, the scattering process, represented generically by T. From a source 
point of view, both processes contribute to the vacuum amplitude 
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Fig. 7.7 Two particles produced by the two earlier sources scatter by the processes 
labeled T and the scattered particles are detected by the two later sources. 


W = - 


K*G°K 


K*G°TG°K, 


(7.65) 


where the last term means 

- J K*{l)G°il - l')T(l', 1")G°(1" - 

= - j (7-66) 

The vacuum persistence amplitude is expressed as 

(0+|0_) = e*^ (7.67) 

for the same reason as before. Because we have a well-defined causal situation, 

= J G°il-l')Kil'), h>tr, (7.68) 


we can write 


^ EE ^^p(r,t)(-t)/t:p. (7.69) 

p V ^ ^ p 
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T 


Fig. 7.8 Diagram representing scattering of a particle off a fixed center. Here the 
heavy vertical line represents the fixed center or nucleus, the dotted line represents 
the exchange of a particle directly between the sources, and the solid lines represent 
the production of a particle by the first source, its scattering off the fixed center, by 
the process represented by the black box T, and finally its absorption by the second 
source. 


Similarly, 



(7.70) 


implies 



(7.71) 


The term iW is the only term which describes the process being considered. 
The vacuum amplitude includes 


^(o+|ip)^*(ip+lV_)(V+)^ « ^(-zii:;)(ip+|V_)(-tii:pO. (7.72) 


P 


P 


So, picking out the coefficients of — —iK^^ we infer 



(7.73) 
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SO if we knew T we could compute the scattering amplitude by taking matrix 
elements in this way. 

Let us analyze the scattering process in more detail, by going into processes 
which occur in a definite time. That is, we break up the extended process, 
represented by the black box in Fig. 17.81 by viewing it as a repetition of 
elementary processes represented by a potential V, as sketched in Fig. 17.91 
To describe this, we extend the idea of the field. The vacuum amplitude is 



V- 


Fig. 7.9 Multiple-scattering resolution of the scattering operator T. Here the di¬ 
agram represents repetitions of elementary scattering processes each described by a 
potential V. 

given in terms of 



(7.74) 


where the quantity in square brackets is regarded as the field, 



(7.75) 


where = f G^K and T'ljP may be thought of as an effective source. In this 
way we get a superposition of effects. Alternatively, we may emphasize the 
last scattering act. 
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where, as above, ip is due to an infinite number of elementary scattering acts. 
This decomposition is self-consistent, because we may write 


iP = iP^+ / G^Vip 


= '^° + J G°Vip° + j G°VG°Vip° + j G°VG°VG°Vip + ... .(7.77) 


So if the expansion makes sense 

T = V + VG°V + VG°VG°V + .... (7.78) 

The differential equation corresponding to Eq. (17.761) is 

= + (7.79) 

or 

- T - e) = X. (7.80) 

So we see, indeed, that V is the potential energy. Thus, we indeed have a 
multiple scattering process. 

This also generalizes the concept of the Green’s function. If we write 


iP 



(7.81) 


where the Green’s function satisfies 

G(l,T) = d(l,l'), 

and the action can be written as 

IT = - y K*iP - J KiP* + J ^ -T-V^ip, 

which in value is equal to 


(7.82) 


(7.83) 


W = -i 


K*GK. 


(7.84) 


The Green’s functions include the possibility of bound states. 

More generally, think of 2-particle scattering, represented in Fig. 17.71 
Single-particle exchange is represented by 


W2 = - 


/ 


K*G°K. 


(7.85) 
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But here, the particles do their thing, and the scattering is represented by 

W4 = -^J {K*G°)i{K*G°)2T{12, 1'2'){G°K)i,{G°K)2>. (7.86) 

Again, the scattering amplitude (Ipilpillpj, Ipa/) is given simply in terms 
of T. Moreover, we want to analyze the scattering in terms of elementary 
processes, as sketched in Fig. 17.101 The terms in the vacuum amplitude of 



T 


Fig. 7.10 Scattering process described in terms of multiple scattering. 


interest are contained in which when expanded contains the fol¬ 

lowing terms with four sources: 

{K*G°){K*G°)TiG°K){G°K) 
= -^l K*{1)K*{1')^{1,2), (7.87) 


- +iW 4 = -^ j 


where 

t/^(l, 2) = t/:°(l)t/:°(2) -b i J (7.88) 

This says that the source detects what comes directly, and what comes from 
the scattering process. On the other hand, we can look at this from the last 
scattering point of view, 
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^'(l, 2) = V'°(l)^°(2) + z J G°G°V^il', 2'), (7.89) 

where V is the measure of the single-scattering process. Again one can iterate. 
Look at the corresponding differential equation, 


d 


dti 


d 


dto 


z/>(1, 2) = K{l)K{2) + f(Fz^)(l, 2), (7.90) 


where a matrix notation is used in the last term. This can be written as 


d 


dti 


d 


zt— -T2]-iV 


dt2 


z/-(l,2) = K(l)iL(2). 


(7.91) 


Now make explicit the time localization of the potential: 

1/(1,2; T, 2') = 5{G - t2)S{ti, - t2,)S{ti - ti,)Vi2{ri - r2, - r2'). (7.92) 

It might also be that the potential is localized in space: 

Vi 2 (ri - r 2 ,ri/ - r 2 ') = (5((ri - r 2 ) - (rn - r 2 '))l^(ri - ^ 2 ), (7.93) 

but we won’t assume this in the following. Then the Green’s function defined 

by 

^'(l, 2) = J G(l, 2 ; T, 2')K{r)K{2), (7.94) 


satisfies 






G = (<5(1 - l')<5(2 - 2')) 


sym : 


(7.95) 


where the subscript denotes symmetrization. Now because of Eq. (17.891) , the 
integral equation for the Green’s function is (symmetrization suppressed) 


G = G°G° -h iG'^G^’VG. 


'0 /^Ot 


(7.96) 


The assumption of an instantaneous V means that we can look at the equal¬ 
time Green’s function, 


G(ri,r2,f;r(,r2,t')i t = h = t2, t' = t[=t2. 
Appearing in Eq. (17.961) is 


(7.97) 


iG°{ri,t;r[,t')G\r2,t-,r2,t') = 

X exp[z(pi • (ri - r() -b z(p2 • (r2 - r2) - 


f (dpi) (dp2) 

J (27r)3 (27r)3 

i{Ei + E2){t - 


t')],(7.98) 


a Green’s function with energy Ei + E 2 . This obeys the equation 
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- Ti - T 2 ^ (tG°G°) = 5(t - t')Siri - r[)Sir 2 - r^), (7.99) 

the conventional Green’s function equation. Then 

{'Ft ~ ~ + Vi 2 {iG), (7.100) 

which is what we would expect. 

Now we write the vacuum amplitude as 

(0+|0_)^ = 1 - i y K*G°K -^J K*K*i;{12) + ... 

^ ^ ^ (7.101) 

where in the last three-particle and higher interactions have been omitted. 
Here we include in the second term only the effects of interaction, through 


x(l,2) = t/;(l,2)-V'(l)V’(2). 


(7.102) 


The integral equation (17.961) written in terms of the last interaction can 
also be written in term of the first interaction, 


G = G°G° -f GiVG°G°. 
Symbolically we can solve for G, 

1 


G = 


1 - G°G°iV 


_ G^G^ 


1 


1 - iVG°G° ’ 


(7.103) 


(7.104) 


which are formally identical. This assumes the instantaneous interaction given 
by Eq. (I7.92|) . Let us write 

G°G° = -iG?+2, iG = G 1 + 2 ; (7.105) 

then our integral equation reads 

Gi+2 = G5^2 + G'5^2^12Gi+2. (7.106) 

which satisfies the differential equation (I7.100L or 

{'Ft - - ^2 - ^ 12 ) Gi+2 = S{t - t')S{ri - ry) 6 {r 2 - r 2 '). (7.107) 

We might be interested in non-instantaneous initial and final states, even 
though we are assuming an instantaneous interaction. Then we could write 


G = iG°(t> - t<)Gi2(t< - <'>)G°(t'> - t'<), 


(7.108) 










no 


7 Nonrelativistic Source Theory 


where 


G^{t — t') ^—i6{r — r') as t —>■ + 0 


(7.109) 


hence the factor of i is supplied by comparison with the instantaneous limit, 
h = t2, t'l = t'2. 


7.2 Bound States 

Sources must be able to create composite structures. We could have started 
with composite particles, after all, so this is a aspect of self-consistency. Write 
the integral equation for G symmetrically. 


G = G°G° + G°G°iV[G°G° + GiVG°G°] = G°G° + G°G°iVG°G° 


+ G°G°iVGiVG°G°. 


(7.110) 


We are still working with the instantaneous interaction approximation. So G 
has the property of depending on two times. Write again 


(7.111) 


G — —iGi+2, V —>■ Vi2. 


The description of the bound states is contained in G 1 + 2 . The part of the 
vacuum amplitude involving G 1+2 is 


(0+|0_) = e-5/^f*'f^*G“GWi2G;+2Vi2G«G«iCif 


(7.112) 


This involves all interactions. 

Now introduce the coordinates for the two-particle system, where we as¬ 
sume that the two particles have the same mass. 



(7.113) 


The Green’s function equation (17.1071) becomes 



r')5(R-R'), 

(7.114) 


where we see the appearance of the total mass M, the reduced mass fi, the 
total momentum P and the relative momentum p. We are interested in the 
motion of the center of mass. Let the internal motion be described by an 
eigenfunction governed by a Schrodinger equation, 
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SO multiplying Eq. (17.1141) by (r) and integrating over r we obtain 

= ,5(t-t')'^(R'-RV*(i’')- (7.116) 

The Green’s function can be expanded, therefore, in terms of a single-particle 
Green’s functions depending on the state n: 

G(r, R, t-r', R', t') = ^ <^„(r)G„(R - R', t - t'X(r')- (7.117) 

n 

Now from Eq. (|7.112l) we identify the effective source for an atomic state 

Kn{'R,t) = J{dr)(j)*{r)Vi 2 {r)G°(R+ ^ - ri,t - ti) 

xG°(R - I - r2,t - t2)K{ri,ti)K{r2,t2)dl d2.{7.118) 

We must be very explicit that En < 0, meaning that the sources are extended: 
They must put out less energy than two free particles would have. So G° does 
not refer to the propagation of a free particle, but rather it propagates an 
excitation which does not get very far. As a result, the Green’s functions 
become real in effect, by virtue of the integration, and 

K*{K, t) = ^j K*K*G°G^Vi2<l)n (7.119) 

is actually the complex conjugate of Eq. (j7.118l) . The consistency of the vac¬ 
uum persistence amplitude 

(0+|0_) = (7.120) 


demands that K* really be the complex conjugate of Kn, 

Perhaps this remark is clarified by looking at the Fourier transform of the 
free Green’s function. 


1 

E-T + ie' 


(7.121) 


When E < 0 there is no singularity and the Green’s function is real. 

The formalism is flexible: it can deal with bound states whether analyzed 
or not, which is especially useful in high energy physics. 

Can we write down a formalism that expresses the dynamics and supplies 
the field equations? The affirmative answer is supplied by writing 


(0+|0_) = e*^, 


(7.122) 


with 
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w = 


-\rrvx-\x*v^i^^ 



(7.123) 


Here we have introduce the 2-particle field x- What does this imply under 
field variations? Varying with respect to ip* gives 


.--r)^ = /w 


■ilj*V{ipip + x), 


(7.124) 


where in the last we see a realistic representation of a single-particle source 
in the interaction with the other fields. Varying with respect to x* yields 




iV 


X = iVipip. 


(7.125) 


This says that x is the part of the two-particle field that has interacted at 
least once, thus 

X = GiVipip. (7.126) 

If we use these equations to evaluate W we recover the previous result (17.1011) . 

We have here two independent field ip and Xi coupled by interaction. This 
is the essence of the many-particle situation in high energy physics. Essen¬ 
tially, we may regard the — ^x*^V'V' ^'S a phenomenological coupling between 
particles. 

This nonrelativistic development of source theory was written down after 
the relativistic formulation was developed, which is explicated in great detail 
in Schwinger’s three volume treatise [Schwinger 1989a'. 











Chapter 8 

Concluding Remarks 


We have traced Schwinger’s development of action formulations from clas¬ 
sical systems of particles and fields, to the description of quantum dynam¬ 
ics through the Quantum Action Principle. In the latter, we have described 
quantum mechanical systems, especially the driven harmonic oscillator. This 
is ahistorical, since Schwinger first developed his quantum dynamical princi¬ 
ple in the context of quantum electrodynamics in the early 1950s, and only 
nearly a decade later applied it to quantum mechanics, which is field theory 
in one dimension—time. At roughly the same time he was thinking about 
quantum statistical systems [Martin 1959) . and it was natural to turn to 
a description of nonequilibrium systems, which was the motivation of the 
time-cycle method, although Schwinger put it in a general, although sim¬ 
plified, context. The time cycle method was immediately applied to quan¬ 
tum field theory by his students, K. T. Mahanthappa and P. M. Bakshi 
[Mahanthappa 1962[ IBakshi 1963] . Then we give a sketch of the application 
of these methods to quantum field theory, based on Schwinger’s 1956 Stanford 
lectures, and to what Schwinger perceived as the successor to field theory. 
Source Theory. The latter appeared shortly after he received his Nobel prize 
in 1963. In fact the present document, which includes components of Julian 
Schwinger’s thinking over nearly four decades, shows that these developments 
proceeding organically, and that the action principle and Green’s function^i] 
played central roles throughout his remarkable career in physics. In fact, it 
has been argued [Mehra 2000] that the first “source theory” paper was in 
fact his most cited one, written in 1951 [Schwinger 1951a] . 
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^ One of Schwinger’s last publications [Schwinger 1993] described the centrality of 
Green’s functions to his life work. 


113 














114 


8 Concluding Remarks 


reported here was used as the basis of lectures in electrodynamics, quantum mechan¬ 
ics, and quantum field theory. 1 especially thank Herb Fried for permission to use his 
transcription of Schwinger’s 1956 lectures as the basis for Chapter 6 here, and Walter 
Becker for his conversion of those notes. 


References 


115 


References 


Bakshi 1963. P. M Bakshi and K. T. Mahanthappa. 1963. “Expectation Value For¬ 
malism in Quantum Field Theory.” J. Math. Phys. 4: 1. 

Dirac 1927. P. A. M. Dirac. 1927. “The Physical Interpretation of the Quantum 
Dynamics.” Proc. Roy. Soc. London A 113 (765): 621-641. 

Dirac 1932. P. A. M. Dirac, V. A. Fock, and B. Podolsky. 1932. “On Quantum Elec¬ 
trodynamics.” Phys. Zeits. Sowjetunion 2: 468. 

Dirac 1933. P. A. M. Dirac. 1933. “The Lagrangian in Quantum Mechanics.” Phys. 
Zeits. Sowjetunion 3: 64. 

Duffin 1938. R.J. Dufhn. 1938. “On The Characteristic Matrices of Covariant Sys¬ 
tems,” Phys. Rev. 54, 1114, doi:10.1103/PhysRev.54.1114. 

Euler 1744. L. Euler. 1744. Methodus Inveniendi Lineas Curvas Maximi Minive Pro- 
prietate Gaudentes. Bousquet, Lausanne and Geneva. 

Euler 1752. L. Euler. 1752. “Investigation of the letter, allegedly written by 
Leibniz,” translated by Wikisource [http://en.wikisource.org/wiki/ Investiga- 
tion_of_the_letter_oLLeibniz]. 

Feynman 1942. R. P. Feynman. 1942. “The Principles of Least Action in Quantum 
Mechanics.” Ph.D Dissertation, Princeton LFniversity, Princeton, NJ. (University 
Microfilms, Ann Arbor, Publications No. 2948). 

Feynman 1949. R. P. Feynman. 1949. “Space-Time Approach to Quantum Electro¬ 
dynamics.” Phys. Rev. 76: 769. 

Feynman 1965. R. P. Feynman and A. R. Hibbs. 1965. Quantum Meehanics and Path 
Integrals. McGraw-Hill, New York. 

Glauber 1963. R. J. Glauber. 1963. “Coherent and Incoherent states of radiation 
field,” Phys. Rev. 131: 2766-2788. 

Hamilton 1834. W. R. Hamilton. 1834. “On a General Method In Dynamics, Part I.” 
Phil. Trans. Roy. Soe. 124: 247-308. 

Hamilton 1835. W. R. Hamilton. 1835. “On a General Method in Dynamics, Part 
H.” Phil. Trans. Roy. Soc. 125:95-144. 

Heisenberg 1929. W. Heisenberg and W. Pauli. 1929. “Zur Quantendynamik der 
Wellenfelder.” Z. Phys. 56: 1-61. 

Keldysh 1964. L. V. Keldysh. 1964. “Diagram Technique for Nonequilibrium Pro¬ 
cesses,” Zh. Eksp. Tear. Fiz. 47: 1515-1527. [English translation: 1965. Soviet 
Physics JETP 20: 1018-1026.] 

Kemmer 1939. N. Kemmer. 1939. “The particle aspect of meson theory,” Proc. Roy. 
Soc. A 173, 91-116. doi:10.1098/rspa.l939.0131. 

Lagrange 1788. J.-L. Lagrange. 1788. Mechanique Analytique. p. 226. 

Mahanthappa 1962. K. T. Mahanthappa. 1962. “Multiple Production of Photons in 
Quantum Electrodynamics.” Phys. Rev. 126: 329. 

Martin 1959. P. G. Martin and J. Schwinger. 1959. “Theory of Many-Particle Sys¬ 
tems.” Phys. Rev. 115: 1342. 

Maupertuis 1744. P. L. M. de Maupertuis. 1744, “Accord de differentes lois de la 
nature qui avaient jusqu’ici paru incompatibles.” Mem. As. Sc. Paris p. 417. 

Maupertuis 1746. P. L. M. de Maupertuis. 1746. “Le lois de mouvement et du repos, 
deduites d’un principe de metaphysique.” Mem. Ac. Berlin, p. 267. 

Mehra 2000. J. Mehra and K. A. Milton. 2000. Climbing the Mountain: The Scien¬ 
tific Biography of Julian Schwinger. Oxford University Press, Oxford. 

Milton 2006. K. A. Milton and J. Schwinger. 2006. Electromagnetic Radiation: Vari¬ 
ational Methods, Waveguides and Accelerators. Springer-Verlag, Berlin. 

Milton 2007. K. A. Milton. 2007. “In Appreciation Julian Schwinger: From Nu¬ 
clear Physics and Quantum Electrodynamics to Source Theory and Beyond.” 
[arXiv:physics/0610054] Physics in Perspective 9: 70-114. 


116 


8 Concluding Remarks 


Noether 1918. E. Noether. 1918. “Invariante Variationsprobleme.” Nachr. Konig. 

Gesellsch. Wiss. Gottingen, Math-phys. Klasse 1918: 235-257. 

Petiau 1936. G. Petiau. 1936. University of Paris thesis, Acad. Roy. de Belg., A. Sci. 
Mem. Collect. 16, No. 2, 1. 

Schrodinger 1926. E. Schrodinger. 1926. “Der stetige Ubergang von der Mikro- zur 
Makromechanik.” Naturwissenschaften 14: 664-666. 

Schwinger 1951a. J. Schwinger. 1951. “On gauge invariance and vacuum polariza¬ 
tion,” Phys. Rev. 82: 664. 

Schwinger 1951b. J. Schwinger. 1951. “The Theory of Quantized Fields. 1” Phys. 
Rev. 82: 914. 

Schwinger 1953. J. Schwinger. 1953. “Theory of Quantized Fields. 111.” Phys. Rev. 
91: 728-740. 

Schwinger 1956. J. Schwinger. 1956. Differential Equations of Quantum Field The¬ 
ory, lectures transcribed by R. Blankenbecler, H. Fried, and J. Tiemann, issued 
by Stanford Research Institute, and transcribed into by W. Becker. 

Schwinger 1960a. J. Schwinger. 1960. “The Geometry of Quantum States.” Proc. 
Natl. Acad. Sci. USA 46: 257. 

Schwinger 1960b. J. Schwinger. 1960. “Unitary Operator Bases.” Proc. Natl. Acad. 
Sci. USA 46: 570. 

Schwinger 1960c. J. Schwinger. 1960. “Unitary Transformations and the Action Prin¬ 
ciple.” Proc. Natl. Acad. Sci. USA 46: 883. 

Schwinger 1960d. J. Schwinger. 1960. “The Special Canonical Group.” Proc. Natl. 
Acad. Sci. USA 46: 1401. 

Schwinger 1961. J. Schwinger. 1961. “Brownian Motion of a Quantum Oscillator.” 
J. Math. Phys. 2: 407. 

Schwinger 1970a. J. Schwinger. 1970. Quantum Kinematics and Dynamics. Ben¬ 
jamin, New York. 

Schwinger 1970b. J. Schwinger. 1970. Particles, Sources, and Fields. Addison- 
Wesley, Reading, MA. 

Schwinger 1973. J. Schwinger. 1973. “A Report on Quantum Electrodynamics,” in 
J. Mehra, The Physicist’s Conception of Nature”. Reidel, Dordrecht. 

Schwinger 1973. J. Schwinger. 1973. Particles, Sources, and Fields, Vol. II. Addison- 
Wesley, Reading, MA. 

Schwinger 1989a. J. Schwinger. 1970, 1973, 1989. Particles, Sources, and Fields, 
Vol. I-III. Addison-Wesley, Redwood City, CA. 

Schwinger 1989b. J. Schwinger. 1989. “A Path to Quantum Electrodynamics” 
Physics Today, February. [Reprinted in Most of the Good Stuff: Memories of 
Richard Feynman. 1993. (Eds. L. M. Brown and J. S. Rigden) AIP, New York.] 
Schwinger 1993. J. Schwinger. 1993. “The Greening of quantum field theory: George 
and I,” arXiv:hep-ph/9310283. Published in Julian Schwinger: The Physicist, 
the Teacher, and the Man, ed. Y. J. Ng. World Scientific, Singapore, 1996, pp. 
13-27. 

Schwinger 1998. J. Schwinger, L. L. DeRaad, Jr., K. A. Milton, and W.-y. Tsai. 1998. 

Classical Electrodynamics. Perseus/Westview, New York. 

Schwinger 2001. J. Schwinger. 2001. Quantum Mechanics: Symbolism of Atomic 
Measurements. Springer, Berlin. 

Sommerfeld 1964. A. Sommerfeld. 1964. Mechanics—Lectures on Theoretical 
Physics, Volume I. Academic Press, New York. [Translated from the fourth Ger¬ 
man edition by Martin O. Stern.] 

Toms 2007. David Toms. 2007, 2012. The Schwinger Action Principle and Effective 
Action. Cambridge Monographs on Mathematical Physics. 

Weyl 1919. H. Weyl. 1919. “Eine neue Erweiterung der Retativitatstheorie.” Ann. 
der Phys. 59: 101-133. 


